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Abstract —Positive semi-definite expressions for the time-domain macroscopic energy density in 
passive, spatially nondispersive, dipolar continua are derived from the underlying microscopic Maxwell 
equations satisfied by classical models of discrete bound dipolar molecules or inclusions of the material 
or metamaterial continua. The microscopic derivation reveals two distinct positive semi-definite 
macroscopic energy expressions, one that applies to diamagnetic continua (induced magnetic dipole 
moments) and another that applies to paramagnetic continua (alignment of “permanent” magnetic 
dipole moments), which includes ferro(i)magnetic and antiferromagnetic materials. The diamagnetic 
dipoles are “unconditionally passive” in that their Amperian (circulating electric current) magnetic 
dipole moments are zero in the absence of applied fields. The analysis of paramagnetic continua, 
whose magnetization is caused by the alignment of randomly oriented “permanent” Amperian magnetic 
dipole moments that dominate any induced diamagnetic magnetization, is greatly simplified by first 
proving that the microscopic power equations for rotating “permanent” Amperian magnetic dipoles 
(which are shown to not satisfy unconditional passivity) reduce effectively to the same power equations 
obeyed by rotating unconditionally passive magnetic-charge magnetic dipoles. The difference between 
the macroscopic paramagnetic and diamagnetic energy expressions is equal to a “hidden energy” that 
parallels the hidden momentum often attributed to Amperian magnetic dipoles. The microscopic 
derivation reveals that this hidden energy is drawn from the reservoir of inductive energy in the initial 
paramagnetic microscopic Amperian magnetic dipole moments. The macroscopic, positive semi-definite, 
time-domain energy expressions are applied to lossless bianisotropic media to determine the inequalities 
obeyed by the frequency-domain bianisotropic constitutive parameters. Subtleties associated with the 
causality as well as the group and energy-transport velocities for diamagnetic media are discussed in 
view of the diamagnetic inequalities. 


1. INTRODUCTION 


As a way of introducing the motivation for and purpose of this paper, begin with the familiar Maxwell 
macroscopic equations that hold in dipolar continua (without current J) to obtain the macroscopic 
Poynting theorem a^ 


with 


J n • [E(r, t) X H(r, t)](i5 = J 


(9D(r,t) 

di 


E(r, t) -|- 


5B(r, t) 


H(r,t) 
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( 2 ) 


D = eoE-FP, B = ;Uo(H-bM) 


( 3 ) 


^ In reading Poynting’s original 1884 paper, one finds that he wrote his theorem in the form (using modern notation and SI units 
with displacement current written as dD/dt instead of edE/dt and written as dB/dt) 



dD ^ dB „ ^ ^ 

'df ' ~dt ' + ■ ^ 





( 1 ) 


where Ev = E + v x B and v is the velocity of the “matter” ^ eq. (7)]. This is not too surprising since Poynting references Maxwell’s 
Treatise in which Maxwell wrote some of his final equations in terms of E + v x B [2], m- The V terms in Q cancel to leave the 
familiar form of “Poynting’s theorem”. 
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where (E, D) and (B, H) are the macroscopic electric and magnetic helds, respectively, P and M are the 
macroscopic electric and magnetic dipolarization densities, respectively, and eo and are the free-space 
permittivity and permeability, respectively. (The term “macroscopic” is used throughout to designate 
electromagnetic fields and sources averaged over macroscopic volumes that contain a large number 
of discrete dipolar molecules or inclusions while having dimensions much smaller than the minimum 
temporal and spatial wavelengths in free-space and in the continuum. An electromagnetic “dipolar 
continuum” or just “continuum” is used throughout to refer to a medium in which fields and sources 
obey the traditional Maxwell dipolarization equations; see (P5|) or (HUIl . We will later find it useful 
to distinguish between “ideal” and “macroscopic” continua.) The unit normal n to the closed surface 
S points out of its volume V so that the left-hand side (and thus the right-hand side) of ([2]) is (as 
will be proven later) equal to the instantaneous electromagnetic power flow entering the volume V ; see 
Fig. m Integrating ([2]) from time to when the maeroscopic fields are zero to the present time t gives a 
macroscopic electromagnetic energy density on the right-hand side of ([2]) equal to 

t 


io 

If it is assumed that the macroscopic electromagnetic energy in a passive polarized continuum is at 
least as great as the energy of the same E and H fields produced in a vacuum, then this vacuum energy 
density is (eo|Ep -|- /io|Hp)/2 and (jl]) yields the inequality 




dt'. 


( 4 ) 
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or simply 
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( 6 ) 


If, on the other hand, one takes the view that all known magnetization is produced by Amperian 
magnetic dipoles (circulating electric current) and thus E and B are the primary fields, then the vacuum 
energy density is (eo|Ep -|- |Bp//xo)/2 and (jH) yields instead of ([5]) and ([6]) 
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to 


9D(r,F) , dB{r,t') 
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( 7 ) 


( 8 ) 


Although it is often argued that ([U]) is the correct inequality for passive continua iHS], m p- 81], 
in fact, neither one of the inequalities in ([6|) or ([8|) are universally valid because the simple hypothetical 
example of a passive material with constitutive relations D = eE and B = /iH with real constant 
permittivity e and permeability /i over the operational bandwidth reveals from ([6|) that 

e > eo , M > |Uo (9) 

and from ([5]) that 

e > eo , 0< fo- (10) 

This would suggest that Q and dSI) may be valid for paramagnetic and diamagnetic materials, 
respectively. Indeed, the main purpose of this paper is to show that, under quite general sufficient 
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Figure 1. Volume V with surface S in dipolar continua with macroscopic polarization densities P and 

M. 

conditions and without assuming any particular constitutive parameters, dH) and ([8|) apply to 
paramagnetic and diamagnetic dipolar continua, respectively. 

To prove these positive semi-definite time-domain inequalities in Q and ([8]), the macroscopic 
Maxwell equations and Poynting theorem are systematically related to the more fundamental 
microscopic Maxwell equations and Poynting theorem from which their macroscopic counterparts are 
derived. Diamagnetic material is defined simply as material with Amperian magnetization produced 
by molecules or inclusions that have no permanent magnetic dipole moments, only induced magnetic 
dipole moments. Paramagnetic material, a term that includes ferro(i)magnetic and antiferromagnetic 
material [H chs. 11 and 12], is defined simply as material with magnetization produced by the alignment 
of “permanent” (sometimes called “intrinsic”) Amperian magnetic dipole moments of molecules or 
inclusions that dominate any induced diamagnetic magnetization^ Notably, the detailed microscopic 
derivation using classical models of electric and magnetic dipoles reveals that the difference between the 
integrals of the power densities in Q and (l8|), what can be called the macroscopic hidden energy, is drawn 
from the reservoir of inductive energy in the permanent (initial paramagnetic) microscopic Amperian 
magnetic dipoles. (The term “permanent” does not imply that the value of the initial current and 
magnitude of the magnetic dipole moment cannot change slightly as the magnetic dipole rotates in an 
external field.) 

Besides the Introduction and Conclusion, the paper contains six other main sections: 

In Section El the mathematically rigorous vector-scalar potential solution is given to Maxwell’s 
differential equations for microscopic electric current and charge density. From these equations, it is 
established that the microscopic Poynting vector integrated over a closed surface in free space equals the 
total instantaneous power crossing that free-space surface. It is further shown that, in passive materials, 
the energy supplied by the fields to the microscopic charge carriers from the time the external fields are 
first applied is always positive semi-definite (nonnegative). 

In Section [3l Maxwell’s continuum equations are derived from the microscopic equations of Section 
El for an ideal continuum, defined by infinitesimal subvolumes of electric and magnetic dipolarization 
that are continuous, and a straightforward proof is given that these same Maxwell equations describe 
the macroscopic dipolarization obtained from averaging the sources and fields of classical models of 
discrete electric and magnetic dipoles. It is shown that a requirement for the equivalence of the ideal 
and macroscopic continua equations is that the surfaces of the macroscopic averaging volumes lie in free 
space and not cut through the discrete dipoles. 

In Section [H the conditions are found (namely that the continuum is effectively spatially 
nondispersive) for the macroscopic continuum Poynting vector integrated over a closed surface to equal 


^ These fundamental definitions of diamagnetism and paramagnetism differ from the more common definitions as materials (or 
metamaterials) with negative and positive susceptibilities, respectively [8]. The more fundamental definitions allow for the possibility 
of both diamagnetic and paramagnetic susceptibilities to be positive or negative at higher frequencies where resonances may occur. 
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the instantaneous power flow across that closed surface. 

Section [5] begins with a derivation of expressions for macroscopic continua energy densities in 
terms of the microscopic energy densities. It is determined that, under sufficient conditions that are 
designated by the term “unconditionally passive”, one of these energy-density expressions remains 
positive semi-definite for diamagnetic media, and a second different energy-density expression remains 
positive semi-definite for paramagnetic media. 

In Section El the important positive semi-definite energy-density expressions for macroscopic 
continua are summarized for paramagnetic and diamagnetic continua and it is shown that the difference 
between these two energy densities equals a “hidden energy” extracted from the energy in the 
“permanent” Amperian microscopic magnetic dipole moments that align to produce the paramagnetic 
dipolarization. 

In Section [71 the positive semi-definite time-domain inequalities of Section [6] are applied to 
paramagnetic and diamagnetic bianisotropic continua to derive inequalities satisfied by lossless 
frequency-domain bianisotropic constitutive parameters. Lastly, subtleties associated with diamagnetic 
continua are discussed, namely, the causality of diamagnetic permeability as well as the group and 
energy-transport velocities in lossless diamagnetic media. 

We employ a classical analysis of macroscopic dipolar continua in which no attempt is made to 
determine the detailed quantum nature of the polarizations and fields of atoms and molecules. It is 
assumed, as in most texts on the electromagnetic fields and equations of polarized media, that since 
most materials below optical frequencies and dipolar metamaterials well below bandgap frequencies 
are adequately described by the classical Maxwell macroscopic equations for a dipolar continuum, the 
microscopic dipoles producing the macroscopic dipolarization can be adequately modeled pragmatically 
by classical electric-charge electric dipoles and Amperian-current magnetic dipoles, irrespective of their 
actual quantum origin. 

There are several reasons for determining nonnegative macroscopic energies from the microscopic 
Maxwell equations for classical diamagnetic and paramagnetic dipolar continua. First, it is one of the 
remaining unresolved problems in classical electromagnetic theory even though the problem is relatively 
easy to state: given a volume of a macroscopic continuum satisfying Maxwell’s dipolar equations and 
illuminated by external fields, are there macroscopic polarization energy densities that never become 
negative for all times after an initial time when the macroscopic fields and polarizations are zero. 
Nonnegative energy expressions provide a means for determining realistic physical limitations, such 
as the upper bounds on the bandwidth and gain of antennas or the lower bounds on antenna quality 
factors. Valuable inequalities satisfied by bulk constitutive parameters as well as by the group and energy 
transport velocities in materials and metamaterials can be obtained from nonnegative macroscopic 
energy expressions. Also, knowing the sufficient conditions for the validity of the nonnegative 
macroscopic energy expressions lends insight into the development and utilization of materials and 
metamaterials that may not be subject to the restrictions imposed by these expressions. 

2. MAXWELL’S EQUATIONS FOR MICROSCOPIC ELECTRIC CHARGE AND 
CURRENT 

We assume that we are dealing with macroscopic continua whose molecules or inclusions can be modeled 
by classical microscopic electric charge and current whose fields can be adequately described by the 
following Maxwell differential equations in SI (mksA) units 


^ , . dh(r,t) 

Vxe{r,t)+ =0 

(11a) 

1 ^ , / N de(r,t) , 

Vxb(r,t) eo =j(r,i) 

/xq at 

(11b) 

V • b(r, t) = 0 

(11c) 

eoV • e{r,t) = Q{r,t) 

(lid) 


where e(r, t) and b(r, t) are the primary microscopic electric and magnetic fields at the position r and 
time t, Q{r, t) and j(r, t) are the microscopic electric charge and current densities, and eo and /xq are the 
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free-space permittivity and permeability, respectively. The primary electric and magnetic fields, e(r, t) 
and b(r, t), are defined by the instantaneous Lorentz force, F(r, t) = g[e(r, t) + v x b(r, t)], that they 
exert on a point test charge q moving with velocity v at the time t and position r. (The point test 
charge q is assumed small enough that it does not disturb the sources of the electric and magnetic fields 
being measured. The sources of e(r, t) and b(r, t) are external to q, that is, they do not include the 
fields at (r, t) produced by the moving test charge q. Also, any radiation reaction force on q is assumed 
negligible, as it will be if the acceleration of the charge is negligible or if the charge q is sufficiently small 
since the radiation force is proportional to the square of q [S].) Note that since there are no microscopic 
polarization densities (p and m) in (llip . it follows that the microscopic electric displacement vector d 
is given by d = eoe and the microscopic secondary magnetic field h is given by h = b//xol3 

The rigorous vector-scalar potential solution to Maxwell’s equations m in the frequency domain 


(with e time dependence) is given by [T2l sec. 2.3.6] 

ba;(r) = V X a^(r) (12a) 

e^(r) = ia;a^(r) - V?/;^(r) (12b) 

with the vector and scalar potentials determined from the expressions 

a^(r)=/xolim [ j^(r')G(|r - r'|)dl/' (13a) 

<5^0 J 

V-Vs 

V’a;(r) = —lim [ Q^{r')G{\r-r'\)dV' (13b) 

eo <5^0 J 

V-Vs 

and the scalar Green’s function given by 

2 ^0 1I 

G(|r - r'l) =kQ = u)y/JI^. (13c) 

47r|r — r I 


The limits before the integrals of (|13p are rigorously required to isolate the singularity in the Green’s 
function at r' = r from the volume integration by a “principal volume” Vs (with maximum dimension 
6) to allow rigorous differentiation of the potentials to obtain the fields in (fT2p [T^ sec. 2.3.6], [T3] . 


2.1. Poynting’s theorem for microscopic electric charge and current 


Let the microscopic charge and current be confined to a region of finite extent so that they can be 
enclosed by a volume V whose surface S lies in free space. Then the instantaneous power Pje{t) 
supplied by the fields to the charge-current within V is simply 

Pje(t) = j j(r, t) ■ e(r, t)dV . (14) 


This is proven from the Lorentz force on a volume charge density moving with velocity v, namely 
dF = £»v(e -I- V X h)dV with dP = v ■ dF = g^v ■ edV = j • edFo (Thus, the magnetic field does not 
imnart enerev to the charee carriers.! With the heln of (lllal) and (lllbl). the Dower in (d can be recast 
in the form of the microscopic Poynting theorem 



■edV 


V 


1 d 


V • (e X h/iio) + (eo|ep |b|Vho) 


dV 


(15) 


^ For metamaterial inclusions with microscopic polarization densities p and m, the current j in Gill would include Op/Ot+V x m, 
and the charge g would include —V • p. Thus, the main results of this paper are applicable to these metamaterials as well within the 
bandwidths that they behave as dipolar continua mm , provided either paramagnetism or diamagnetism dominates the macroscopic 
magnetization produced by the inclusions. For example, metal-dielectric inclusions in metamaterial-array continua would be described 
by the electric-dipole/diamagnetic energy expressions. 

^ The Lorentz force on gdV is exerted by the fields of all sources except the fields of gdV itself. However, the self-fields of an 
infinitesimal volume element of continuous volume density of charge and current approach zero as dV —>■ 0 and thus e and b in this 
volume-element Lorentz force are the total electric and magnetic fields in Maxwell’s equations GD- 
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or by means of the divergence theorem as 

J n • (e X h/fj.o)dS = - J edV - J (eo|ep + |bp//xo) dV (16) 

S V V 

with n denoting the unit normal to S pointing out of the volume V. 

As any one of an infinite number of test examples to which (jl6p can be applied, suppose that 
the volume V contains a charged capacitor C that is shorted at t = 0 through an inductor L. Even 
if the capacitor, the inductor, and the wires are lossless conductors, this LC circuit will not oscillate 
forever because it will radiate. At the time t = T, assume that the radiation has practically reduced the 
oscillations of the RL circuit to zero and that the surface S is chosen far enough away that the radiated 
fields have not reached S at the time T. Then (jl6|] becomes after integrating over time from 0“ to T 

T 

\ J eo|e(r,0")pdR = y J 3 {r,t) ■ e{r,t)dVdt + ^ J [eo\e{r, t)\‘^ + \h{r, t)\‘^/ fio] dV (17) 

V 0- V V 

where 0“ denotes the time just before t = 0. The left-hand side of (jl7p is equal to the electrostatic 
energy stored initially by the capacitor, a result that is rigorously proven from the electrostatic equations 
under a quasi-static charging of the capacitor plates. As explained above, the first integral on the right- 
hand side of (fTTP is the energy supplied by the fields to the charge carriers of the current. (It is zero 
for a stationary lossless conductor and greater than zero for a lossy conductor.) Therefore, the second 
integral on the right-hand side of (jl7J) is an additional energy required by conservation of energy. In 
other words, /y[eo|e(r,t)p -|- \h{r,t)f//Jo]dV/2 is the free-space energy stored in the electromagnetic 
fields in the volume V not only for statics but also for time dependent fields at each instant of time t. 
With this established, we see that the right hand side of (fTH]) . and thus the left-hand side of (fTH]) . is the 
time rate of change of the total energy leaving the volume V. In particular, the integral of the time- 
dependent microscopic Poynting vector, that is, fi- [e(r, t) x b(r,t)/^o]) over a closed surface S lying 
in free space is the total instantaneous electromagnetic power leaving that closed surface. Alternatively, 
the total instantaneous electromagnetic power P{t) entering a closed surface S that lies in free space is 
given by 

P{t) = - J n - [e(r, t) x b(r, t)/yLo]dS 
s 

= j j(r, t) ■ e(r, t)dV + J [eo|e(r, t)\‘^ + |b(r, t)| VH dV . (18) 

V V 

This result is true whether or not the microscopic charge-current inside S forms polarized material, 
provided S lies in free space outside the material. However, this result does not imply that n- (e x b//Jo) 
is necessarily the instantaneous free-space power flow per unit area in the direction of fi. In other words, 
the integral of n-(exb//io) over an open free-space surface does not necessarily equal the electromagnetic 
power that flows across that open surface. 

2.1.1. Energy supplied to the charge carriers 

The instantaneous power supplied by the electromagnetic fields to the charge carriers (usually electrons 
and positively charged nuclei) producing the microscopic charge and current (^),j) that generate the 
fields e and b under consideration in V is given in (jl4p . Expressing Pjeit) as the time derivative of an 
energy Wje{t), that is, Pjeit) = dWje{t)/dt, and integrating the power in ([Till from an initial time to in 
the past to the present time t, we obtain 

t 

Wjeit)-Wjeito) = J j 3{r,t') ■ e{r,t')dVdt'. (19) 

to V 
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This is the energy supplied to (work done on) the charge carriers of the current j in y by the electric 
field e produced by the charge carriers during the time interval t — to. Assume that the applied 
electromagnetic fields are zero until after the time to so that the electric field does no work on the 
charge carriers until after to. If the charge-current of the dipoles is produced by bound charge carriers 
(charge-current comprising electric and magnetic dipoles with negligible translational displacement and 
negligible changes in their number per unit volume), then this work, Wje{t) — Wje{to), done on the 
bound charge carriers equals the change in kinetic energy of the charge carriers supplied by e plus 
any change in kinetic energy of the charge carriers caused by forces other than the e-field force on the 
charge carriers^ These other forces are fundamentally electromagnetic forces from the fields of the host 
material, namely from the fields of the atomic and molecular charges other than the charge carriers 
that produce e and b. They are the forces, for example, involved in the conversion of the kinetic energy 
of the charge carriers into heat (out-of-band energy), or into the kinetic-potential-heat energy of the 
springs in a compressed-spring model of electric dipoles, or into the kinetic-potential-heat energy of 
expandable conducting wires in a wire-loop model of Amperian magnetic dipoles. 

If the macroscopic continuum is passive in that (i) there is negligible initial kinetic energy of 
the bound charge carriers within the operational bandwidth (that is, negligible in-band initial kinetic 
energy) that can be decreased, (ii) there is negligible in-band net energy transfer from the host material 
to the charge carriers (that is, in-band energy cannot travel through the host material, except by means 
of e and b, from one region of the continuum to the other), and (iii) there are no auxiliary sources (such 
as chemical reactions changing the e and b fields so as to add to the energy of the charge carriers, or 
initial kinetic-potential energy stored in the springs of a compressed-spring model of electric dipoles, 
or in expandable conducting wires of “permanent” Amperian magnetic dipoles) that can release energy 
upon excitation by applied fields to serve as an active source of internal energy increasing the energy of 
the charge carriers, then Wjeft) — Wjefto) > 0. In other words, if there is no initial in-band mechanical 
energy, and no internal active sources of energy, and no nonlocal transfer of in-band mechanical energy, 
then the material is passive such that Wje{t) — Wjefto) > 0. Choosing the arbitrary, constant initial 
energy Wje{to) equal to zero, we have 

t 

Wje{t) = j J i{r,t') ■ e{r,t')dVdt'>0 (20) 

to V 

for bound charge carriers in passive material that can be lossy as well as lossless, linear or nonlinear. 
Equation (I20p can be taken as the mathematical definition of passivity. A more restrictive “unconditional 
passivity” will be defined later in Section [5l 

The current j(r, t) in (j20p can include conduction current (in addition to the current produced by 
the local motion of the bound charges) if the changes in the energy transferred by the drift velocity 
of the conduction charges are negligible compared with the heat energy produced by the conduction 
current (as is normally the case). 

3. MAXWELL’S EQUATIONS FOR DIPOLAR CONTINUA 

In order to explain how the microscopic electric and magnetic fields of discrete electric and magnetic 
dipoles should be averaged to obtain the traditional macroscopic Maxwell equations for dipolar media, 
it is helpful to first determine the ideal continuum Maxwell equations directly without assuming the 
existence of discrete dipoles. This direct continuum approach for deriving Maxwell’s dipolarization 
equations originates with Maxwell himself mm- 

^ Each differential element of charge gdV could radiate energy and thus experience an irreversible self-radiation-reaction energy as well 
as reversible self-radiation-reaction (often called the “Schott acceleration energy”) and self-electromagnetic-momentum energies [^. 
However, since all these self-energies are proportional to the square of the charge, that is, {gdV)^, they are a higher-order differential 
than dV and thus they can be ignored in the volume integral of ()19ll . (This does not imply that the integrated self-force/momentum 
and self-power/energy for a fixed amount of charge moving as a relativistically rigid charged particle is negligible [^.) 
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3.1. Maxwell’s equations for ideal dipolar continua 

The essence of the derivation of the ideal continuum Maxwell equations is to first determine the vector 
and scalar potentials in free space outside the differential volume elements with dipole moments and then 
to mathematically define these potentials within the polarization source regions by simply applying the 
same mathematical expressions within the source regions. These Maxwell fields defined mathematically 
within the polarization source regions are then related to fields that can be measured in free-space 
cavities formed by removing an infinitesimal volume of polarized material without altering the remaining 
polarization. 

Specifically, the vector and scalar potentials in free space outside a distribution of differential 
volume elements of electric and magnetic dipolarization, Pi^dV and are given in the frequency 

domain (with time dependence) by [Ml secs. 9.2-9.3] 

a^(r) = fioJ [-MP^(r')G + M^(r') x V'G] dV = fio J [-iwP^(rO + V' x M^(r')] GdV (21a) 

V y 

^uj(r) = - [ P..(r ) • V'GdV' =-- / V' • P^(r')GdP' (21b) 

eo y eo y 

V V 

where G is given in (|13cl) and the volume-element contributions are integrated over a volume V whose 
surface S lies in free space so that it encloses all possible equivalent electric surface charge — h' • P^^ 
and equivalent electric surface current x n' (surface delta functions in — V' - P^^ and V' x M^^). 
These equations can be rigorously derived from m and the definition of infinitesimal dipole moments 
determined by electric charge separation and circulating electric current (Amperian magnetic dipoles) 
outside the sources (r / r') of these dipole moments, PojdV and 'M.^dV. We can express P^ and 
in terms of bound charge and current densities and as 

^ dV (22a) 

Ay 

M^(r) = ^Ihn^ J r' X j^b(r + r')dV' (22b) 

Ay 

with the position vector r inside AV, whose surface AS does not intersect any of the charge and current 
densities that produce the electric and magnetic dipolarization, and f^y gujb(^)dV = 0. The magnetic 
and electric fields in (fTTI) outside this dipolarization are given in the frequency domain as in (fT^ 

b^(r) = V X a^(r) (23a) 

e<^(r) = ia;a^(r) - VV’a;(r). (23b) 

Next let us simply define mathematically the vector and scalar potentials inside as well as outside 
the sources P^^ and by the same integrals as in (|2ip . namely 


A^(r) = ^0 lim f [-*wPtj(r') -|- V' x Maj(r')] GdV' 

~'^V-Vs 

qfjr) = lim [ V' • Pjr')GdV' 
eo <5^0 y 
v^-y^ 

with the mathematically defined electric and magnetic fields given from (I23p as m 

Baj(r) = V X Aa;(r) 

E^(r) = fa;A^(r) - VT^(r) 


(24a) 

(24b) 


(25a) 

(25b) 
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where we have now used the symbols E^, A^, and rather than and ipuj because B^, 

and defined mathematically by (p^ ~ (f25]l are not necessarily equal to the microscopic fields 
bo;, and in the source regions of P^j and M^. The principal-volume limits are reinstated 

in the integrals of (j24[) because, unlike in the integrals of (j21[) . now r can lie in the source regions of 
P(^ and Mtj and thus the singularity of the Green’s function G(|r — r'|) at r = r' has to be excluded. 
Nonetheless, the principal value limits must be retained in (1240 to get the correct values of the fields (and 
their spatial derivatives) in the source regions in terms of P^^ and by means of the differentiations 
in and subsequent use of integral identities m- 

Comparison of (I24j) - (I250 with (jl20 - (jl3p . which satisfy the frequency-domain version of the Maxwell 
equations in m, reveals that the fields and polarizations of the ideal continuum also satisfy the 
frequency-domain version of the Maxwell equations in (fTTI) but with E^j and B^ replacing and b^, 
respectively, and with —iujPaj + V x and —V • P^j replacing and g^j, respectively. Consequently, 
the ideal continuum dipolarization fields satisfy the Maxwell equations 



(26a) 

, (9E(r, t) 9P(r, t) 

t) « = gj +VxM(r,t) 

(26b) 

V • B(r, t) = 0 

(26c) 

eoV • E(r,t) = -V • P(r,t) 

(26d) 


where the frequency-domain equations have been converted back to the time domain by taking the 
inverse Fourier transform. In other words, the mathematically rigorous solution for all r to the 
frequency-domain version of (j26p is given by (I24n - (j25p . If the surface (S' of F containing a polarized 
chunk of continuum material lies in free space outside the chunk, the integrations in (|24p include the 
equivalent electric surface charge and current (— n' • P^^ and x n') by means of the surface delta 
functions in —V' • P^ and V' x M(^; whereas if S lies just inside the chunk of material, these surface 
contributions are not included. To the right-hand sides of (I26bp and ()26dh can be added free current 
J/(r,t) and free charge pf{r,t), respectively. 

With the secondary fields D and H dehned by the constitutive relations 


D(r, t) = eoE(r, t) + P(r, t) 

(27a) 

H(r,t) = —B(r, t) — M(r, t) 

(27b) 

li-O 

the equations in (I26p can be rewritten as 


VxE(r,t) + ^^^ = 0 

(28a) 

VxH(r,t) -J/(r,t) 

(28b) 

V • B(r, t) = 0 

(28c) 

V-D(r,t) = pf{r,t). 

(28d) 


The E and B fields, which have been defined mathematically in the source regions of electric 
and magnetic polarization P and M can be related to measurable fields inside a cavity formed by 
instantaneously removing an infinitesimal volume of P and M without disturbing the remaining 
polarization. For a circular-cylinder infinitesimal volume of length 2b and radius a, whose axis is aligned 
with P, the electric field in the free space at the center of the cylinder differs from the mathematically 
defined electric field by the amount m 


AE = 


eo 



y/a? + h‘^ 


(29) 
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For an infinitesimally narrow cylinder (a/6 —>■ 0), AE = 0 and the narrow-cylinder (nc) cavity 
electric field equals E, the mathematically defined electric held, that is 


E 


nc 

c 


= E 


(30) 


where E“ is measurable (in principle) in the free space of the cavity and provides a method for 
determining the primary electric held E. 

If 6/o —)• 0 so that the cylinder becomes a disk, AE = P/eo and the thin-disk (td) cavity electric 
held is given by 

E/d = E + — = — . (31) 

eo eo 

Thus measuring E*^ provides a method for determining the secondary electric held D (often called the 
electric displacement because it is produced by charge separation or “displacement”, as in a capacitor). 
Similarly, for an inhnitesimal circular cylinder aligned with M [13] 


AB = -/UqM 


6 

+ 6 ^ 


(32) 


so that 

B)T = B - = /ioH (33) 

Bf = B . (34) 

Thus the cavity magnetic helds provide a way to measure the mathematically dehned primary magnetic 
held B (often called the magnetic induction because the time derivative of its hux through an open 
surface induces an electromotive force around the edge of the surface, as in an inductor) and the related 
secondary magnetic held H. 


3.2. Maxwell’s equations for macroscopic dipolar continua 

The ideal dipolar-continuum Maxwell equations in (|26p or (|28p have been derived assuming that each 
differential volume element of PdE and MdP have the same properties as the continuum as a whole, 
that is, they are simply inhnitesimal chunks of a dipolar continuum with no free space outside the 
surfaces of the chunks except at the outer surface of the continuum. However, most natural materials 
and many metamaterials are comprised of discrete molecules or inclusions separated from one another 
by a hnite distance in free space. If the operational temporal and spatial bandwidths of the externally 
applied signals are low enough that these molecules exhibit only electric and magnetic dipole moments 
with no signihcant higher order multipole moments, and their average helds vary slowly over electrically 
small macroscopic volumes AV containing many molecules or inclusions, we will now show that these 
average helds obey the same Maxwell equations as in (IMp or ([28]). 

To simplify the derivation, let AV be a spherical electrically small volume containing a large number 
of the discrete electric and magnetic dipoles and assume that the surface AS of AV lies in free space 
without cutting through any of the dipoles^ (We can imagine displacing slightly a few of the dipoles 
or part of the spherical surface to realize this assumption; the fractional error in the macroscopic 
polarization densities and helds introduced by this displacement is on the order of the ratio of the 
diameter of the dipoles to the diameter of the spherical macroscopic volume AH.) The average electric 
and magnetic dipole moments of the bound charge and current Q}, and within AH are given by 

P(r,t) = ^ / ^b{r + r',t)r'dV' (35a) 

AV 

® The polarization densities and fields defined by averaging (at each instant of time t) over the macroscopic volumes AV containing 
discrete numbers of dipoles will not be perfectly continuous functions of position r as individual dipoles are included or not within 
AV as the center r of AV changes position ESI pp- 2-3], [It]. These microscopic discontinuities associated with macroscopic- 
volume averaging over discrete numbers of dipoles can be sufficiently smoothed by various mathematical techniques to allow the curl 
and divergence of the fields in Maxwell’s differential equations to be well-defined. Alternatively, the volume definitions of curl and 
divergence m can be applied directly to the macroscopic volumes. 
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M(r,t) = ^^ / ^ 3 b{r+ r',t)dV' (35b) 

AF 

where r is the position of the center of the spherical volume AV. We can use the same continuum 
polarization symbols P and M defined in (1221) (after converting (1221) to the time domain by taking the 
inverse Fourier transform) for the average macroscopic dipole moments in (|35p because it is assumed that 
the diameter of AV in (|35p is much smaller than the smallest significant wavelength in the operational 
temporal and spatial bandwidths and, thus, the macroscopic volume AV can effectively be used in place 
of the infinitesimal differential volume element dV. It is assumed that the total bound charge in each 
AV is zero so that its electric dipole moment is independent of the position of the origin r with respect 
to AV. The magnetic dipole moment of the bound current in AV is dependent upon the position of the 
origin r, unless the electric dipole moment is zero. However, for r located within AV, the dependence 
of M on position is minimal if AV can be made much smaller than the smallest significant wavelength 
and still contain many discrete dipoles. 

It will now be shown that the electric and magnetic fields averaged over the spherical macroscopic 
volume AV are approximately equal to the average fields within AV of the corresponding ideal 
continuum. First, consider the average fields in AV produced by the dipoles contained inside AV. 
For a finite number N of translationally stationar}0 discrete electric and magnetic dipole moments p* 
and rrij inside the spherical macroscopic volume AV (centered at the position r at time t and whose 
surface AS lies in free space not intersecting any of the dipoles), the P and M in (|35p can also be 
expressed as 

1 ^ 

P(r,i) = ^X]P* 

2=1 

I ^ 

M(r, t) = — ^ m*. (36b) 

2 = 1 

Let the microscopic electric field of any one of the electric-charge electric dipoles pi be denoted by 
ej(r,t). Then the average electric field of this dipole inside the spherical volume AV is defined by 

= ^ j ei{r',t)dV'. (37) 

AV 


For operational bandwidths that are not so large that the diameter of the macroscopic volume can be 
a small fraction of a minimum wavelength (and yet contain many dipoles), the electric field Oj over 
the volume AV is dominated by the quasi-electrostatic fields and, thus, the average in (l37|l is simply 
PaTi ~ “Pi/(3eo) [HI sec. 4.1]. Summing over the N electric dipoles and using (|36ap gives 

k: = ( 38 ) 


where E™® is the average over AV of the quasi-electrostatic field of the dipoles inside the AV centered 
at the position r at time t. Similarly, for the B field averaged over the spherical AV produced by the 
Amperian magnetic dipoles rrij inside AV, one finds with the help of (|36bp that [141 sec. 5.6] 


ins _ 2//oM(r, t) 


(39) 


The 2/3 factor in (f39P replaces the —1/3 factor in (f38p because the magnetic dipoles are produced by 
circulating electric currents rather than equal and opposite magnetic charges. 

The average over AV of the quasi-static electric fields produced by the magnetic dipoles inside 
AV is negligible compared to the average of the quasi-electrostatic fields E™® produced by the electric 

^ Translating dipoles can be relativistically transformed at each instant of time to a sum of stationary electric and magnetic dipole 
moments. 
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dipoles inside Al^. Likewise, the average over AV of the quasi-static magnetic fields produced by the 
electric dipoles inside AV is negligible compared to the average of the quasi-magnetostatic fields B™® 
produced by the magnetic dipoles inside AV. 

Second, consider the fields averaged over AV produced by the dipoles outside AV. The average 
over the spherical volume AV of the quasi-electrostatic fields from the electric dipoles outside AV is 
simply equal to the value of the quasi-electrostatic field at the center of the free-space spherical cavity 
AV] that is [HI sec. 4.1] 

= (40) 

where r is at the center of the sphere. Since the center of the sphere is much further away from the 
dipoles outside AV than the average distance between the dipoles, the discrete dipoles outside AV can 
be approximated by an ideal continuum for the sake of determining E°“*(r, t). Moreover, this ideal 
continuum approximation becomes even better for the dipoles so far away from AV that their radiation 
fields become appreciable to or much larger than the quasi-static fields. Then these distant fields will 
not vary significantly over the electrically small volume AV and their average fields will also equal the 
fields at the center of the spherical AV. Consequently, the electric field of the dipoles outside AV 
averaged over AV is equal to a good approximation to the averaged electric field from the continuum 
outside AC; that is 

Eout ^ J.OUt (4;^^ 

where E°'^* denotes the average electric field in AC produced by the ideal continuum lying outside AC. 
Similarly, 

Bout ^ gout ^ ( 42 ) 

Adding the inside-dipole average fields in (|38]) - (|39]1 to the outside-dipole average fields in (f4T]l " (f42]l 
for the spherical AC, yields the average fields 


_ -[7 out 

flv — li/ — 


Fjvt) 

3eo 


(43a) 


„ T^out , 2/ioM(r,t) 

-t>av — xj + 


(43b) 


of the discrete dipolar material for bandwidths small enough that AC can be quasi-static in size and yet 
contain many dipoles@ Now, as explained above, E°'^* and B°“* are equal to the average electric and 
magnetic fields in the spherical cavity AC produced by an ideal continuum lying outside AC. These 
fields are generally called the spherical cavity fields of an ideal continuum and are given in terms of the 
total fields of the ideal continuum as nails] 


which implies from (1431) that 


gout _ g 1 P('') 0 

3eo 

(44a) 

out _ g 2//oM(r, t) 

3 

(44b) 

Eav(r,f) = E(r,f) 

(45a) 

Bav(r,f) = B(r,f). 

(45b) 


In other words, the macroscopically averaged fields in discrete dipolar materials or metamaterials are 
approximately equal to the mathematically defined ideal continuum fields if the macroscopically averaged 
dipolarizations are approximately equal to the continuum dipolarizations and the minimum temporal 
(free-space) and spatial (medium) wavelengths (call them A™™ and A™™, respectively) are much larger 
than the average separation distance (d) between the discrete dipoles. Thus, these macroscopically 

® The P/(3eo) term in Il43all and Il44all also appears in the local-field expression of the Lorentz-Lorenz relation. However, the local 
field is the electric field at one electric dipole produced by all the other electric dipoles in random or highly symmetric lattices such 
as cubic lattices ch. 16]. It is not an average of the dipole fields over a volume and for asymmetric lattices the P/(3eo) term no 
longer determines the local field. 
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averaged fields and sources obey the same Maxwell equations in [^) or (2^ as the mathematically 
defined ideal continuum fields and sources, namely 


VxE(r,i) + ^^^ = 0 (46a) 

V X H(r, t) - = J/(r, t) (46b) 

V • B(r, t) = 0 (46c) 

V • D(r,t) = p/(r,t) (46d) 

where the free macroscopic electric charge and current densities are given in terms of the free microscopic 
discrete electric charges qi and current densities as 


P/(r,i) = (47a) 

2 = 1 

1 ^ 

J/ (r, t) = ^ X] (47b) 

2 = 1 

in which Vj is the velocity of the charge qt. The macroscopic constitutive relations follow from (1270 as 

D(r, t) = eoE(r, t) + P(r, t) (48a) 

H(r, t) = — B(r, t) - M(r, t) . (48b) 

ho 

If we denote the maximum temporal and spatial wave numbers by = 27r/A™™ and = 
27r/A™™, respectively, then the discrete dipolar medium can be considered a macroscopic continuum 
that approximates an ideal continuum if both k^^d <C 1 and k^^d <C 1. This result, which is one of 
the cornerstones of classical electromagnetics, has also been confirmed in recent articles that rigorously 
analyze three-dimensional cubic arrays of inclusions, which are separated from one another in free space, 
taking into account both temporal and spatial dispersion [iniiiiiEni. It is emphasized that this result 
is derived and, in general, holds only if the surfaces of the infinitesimal volumes used to define the 
dipolarization densities in ideal continua and the surfaces of the finite-size macroscopic volumes used 
to define the dipolarization densities in discrete-dipole macroscopic continua he in free space so as not 
to intersect any of the charge and current that produces the electric and magnetic dipoles. 

I have not been able to find elsewhere the relatively simple, yet rigorous proof given here that 
averaging the microscopic fields of discrete dipoles over macroscopic volumes leads to the same Maxwell 
equations and cavity fields obtained for mathematically defined fields in ideal continuous dipolar media. 
For example, the widely referenced “theory of electrons” used by Lorentz [21] and Rosenfeld [22j to 
obtain the macroscopic Maxwell equations begins by taking the spatial average of the microscopic 
equations in (fTT|) over fixed macroscopic volumes V that cut through the charges and dipoles to get 


V X 5(r,t) + = 0 

(49a) 

VxB{r,t) eo -J{r,t) 

ho ot 

(49b) 

V ■B{r,t) = 0 

(49c) 

eoV • S{r,t) = p{Y,t) 

(49d) 

where the macroscopic average of e(r, t), for example, is defined as 


^(r,t) = ^ / e(r+ r',t)dV' 

(50) 

V 
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and similarly for the other fields and source densities. The macroscopic equations in (I49p hold exactly 
for any size volume V (not just electrically small quasi-static volumes) as long as the limits of integration 
of the volume V in (|5n]l stay fixed for each r and thus the surface S of V, unlike the surface AS of the 
electrically small quasi-static macroscopic volumes AV used above, will cut through the charges and 
dipoles as r is varied. Then, for sufficiently small V, both Lorentz and Rosenfeld claim that the average 
“bound” current and charge on the right-hand sides of (I49bp and (I49dp can be approximated by 


and 

where 


*7(r, i) ~ + V X AT(r, t) (51) 

p{r,t) ^-V ■V{r,t) (52) 

If 1 

'^(rT) ~ y J Qb{r+ r',tydV'^-'^Pi (53a) 

V 

\ f 1 

M{r, ~ ^ y r' X jb(r r', t)dV' ~ y X] 

V 


with Pi and riij equal to the electric and magnetic dipole moments of the molecules within V. However, 
if one inserts the integral definitions of 7^(r,t) and A4(r, t) from (j53p into ()51l) and then the resulting 
J' (r, t) into (]49bl) . one finds with the use of vector-dyadic identities that t) = (1/V) /y jfe(r-|-r', t)dV' 
holds in general only if the surface of V lies in free space and does not cut through jfe(r -|- r',t); that 
is, the V in (l53P cannot be the same as the V in (l50P and thus the validity of equations (|3^ - (f53]) is 
uncertain. This is dramatically illustrated by choosing gi, = —V • T* and jj, = d'P/dt -|- V x M. in (l53P . 

If one ignores this uncertainty and chooses the same V in (|53p as in (j49l) - (j50p . the Lorentz-Rosenfeld 
macroscopic equations in (f4^ - (f53P (which are also obtained by Van Vleck [23l sec. 3]) can lead to large 
variations in the average macroscopic polarizations and fields over microscopic distances (not just the 
microscopic variations, mentioned in Footnote 6, associated with AV). Consider, for example, the 
surface S cutting through discrete electric dipoles composed of extremely large equal and opposite 
charges separated by extremely small distances. Then both the electric polarization and electric field 
averaged over the macroscopic volume V can have extremely large variations over distances comparable 
to the charge separation of the dipoles, and these variations approach infinite values for point dipoles 
with nonzero dipole moments. Moreover, with S cutting through the dipoles, the average of the 
microscopic fields within a cavity V will not generally be a good approximation to the average cavity 
fields of the ideal dipolar continuum equations. 

These large unphysical fluctuations in the Lorentz-Rosenfeld macroscopic equations for dipolar 
continua can be reduced by using smooth test functions or ensemble averaging instead of fixed 
macroscopic volumes V that cut through the dipoles [IT], pa chs. 5-6], pa sec. 6.6]. Nevertheless, 
once the smooth macroscopic Maxwell equations are found, it is still necessary for many theoretical, 
numerical, and experimental purposes (such as deriving the macroscopic power and energy relations 
from the microscopic equations) to determine the conditions under which the cavity fields of these smooth 
ideal continuum equations are to a good approximation equal to the average of the microscopic cavity 
fields of the discrete dipoles. This determination entails introducing macroscopic volumes AV, whose 
surfaces AS do not intersect the microscopic dipoles, and performing a proof similar to the one given 
above in the main body of this section. The necessity of using macroscopic averaging volumes with 
surfaces that do not intersect the sources to obtain physically and mathematically robust, unambiguous 
macroscopic Maxwell continuum equations is further confirmed by the rigorous analysis of spatially 
dispersive metamaterial arrays mm- De Groot and Suttorp, in their book on the foundations of 
electrodynamics p5], average over discrete undivided “stable groups” of point charges that contain 
electric and magnetic dipoles (and higher order multipoles). Also, Landau and Lifshitz, in defining 
macroscopic polarization |26l secs. 6 and 29], state that the surfaces of the macroscopic volumes must 
enclose the dipoles but nowhere enter them. 
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4. POYNTING’S THEOREM FOR DIPOLAR CONTINUA 

Having determined Maxwell’s equations for both an ideal and a macroscopic continuum, we shall derive 
Poynting’s theorem in a way that reveals that the integral of the Poynting vector over a closed surface 
equals the instantaneous power flow in these two types of continua. 

4.1. Poynting’s theorem for ideal dipolar continua 

It was proven in Section 12.11 that the total electromagnetic power at time t entering a closed surface S 
lying in free space is given by the integral over S of the microscopic Poynting vector; specifically from 

m 

P{t) = — J h - [e(r, t) X h(r, t)]dS (54) 

s 

where n is the unit normal pointing away from the volume V enclosed by S, and h has been substituted 
for b//io because S is in free space. We showed in Section IXTl that the fields in an ideal dipolar continuum 
satisfied the Maxwell equations in (|28l) . From these equations it follows that the components of E and 
H tangential to an interface between free space and the polarized material are continuous, provided 
there are no equivalent polarization surface currents at the interface of the polarized material, that 
is, no delta functions in P or M at the interface. Such delta functions can occur for material with 
constitutive parameters that have zero or infinite values m- Even if we ignore these extreme values as 
being unrealizable, finite ( 7 ^ 0 or 00 ) constitutive parameters that are strongly spatially dispersive can 
also exhibit delta functions in P and M at the interface |10llllj . (The usual Poynting vector in strongly 
spatially dispersive material does not generally represent energy flow [26l p. 361], [28].) If we assume 
low enough frequencies that spatial dispersion is negligible, then the tangential components of E and 
H are continuous across the free-space/dipolar-material interface. This implies that if we remove an 
inhnitesimally thin shell of material containing S (everywhere that the closed surface S passes through 
a dipolar continuum), the value of n • (E x H) will be continuous across the free-space shell. Moreover, 
since the continuum fields equal the microscopic fields in the resulting free-space shell, that is 

n • (E X H) = n • (e X h) (55) 

in the free-space shell, one finds from (|54p that the total instantaneous power flowing across a closed 
surface S in an ideal dipolar continuum is given by the integration of the continuum Poynting vector 
P{t) = — f n ■ [E(r, t) X H(r, t)]dS', or with the aid of equations (|28al) and (|28bp 
s 

P(t) = — J h- [E(r, t) X H(r, t)]dS = J 

s V 

Thus, we have determined that in an ideal continuum the integral of the continuum Poynting vector 
over a closed surface S' of a volume V exactly equals the instantaneous power flow across that surface, 
provided the spatial dispersion is not strong enough over the operational bandwidths to produce delta 
functions in P and M at a hypothetical interface S between the polarized material and free space. 
Furthermore, the surface integral of the continuum Poynting vector equals the corresponding volume 
integral of the helds on the right-hand side of (|56p . If the material is linear and characterized by a 
frequency independent permittivity e(r) and permeability ^(r) over the operational bandwidth, then 
(1561) becomes 

P{t) = - j n • [E(r,t) X H(r, t)]fiS’= J [e(r)|E(r, t)|2-h ^(r)|H(r, t)p] dP (57) 

S V 

which implies that in a nondispersive, scalar, linear magnetodielectric material, [e(r)|E(r, t)p -|- 
^(r)|H(r, t)p]/2 is the electromagnetic energy density stored in the E and H fields. 




dt 


dt 


dV. 


(56) 
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It should be noted that in some contexts bianisotropic media are considered to be spatially 
dispersive because the electric and magnetic fields in bianisotropic media are proportional to the curls 
of the magnetic and electric fields, respectively. However, since this weak spatial dispersion maintains 
continuous n • (E x H) across interfaces between the bianisotropic continua and free space, for the 
purposes of this paper, we can include the weak spatial dispersion of bianisotropic media within the 
context of spatially nondispersive media; see Footnote 16. 


4.2. Poynting’s theorem for macroscopic dipolar continua 

The argument leading to Poynting’s theorem in an ideal dipolar continuum and its interpretation in 
terms of power flow can be applied with a couple of modifications to a macroscopic dipolar continuum 
comprised of discrete electric and magnetic dipoles. The first modification is that the volume V with 
surface S to which Poynting’s theorem is applied should be no smaller than a macroscopic volume AV 
which is electrically small but still contains so many dipoles that the averaging of the fields over AV 
gives to a good approximation the continuum fields satisfying Maxwell’s equations in (12611 or ( 12811 and 
ii; see Section Ea 

The second modification involves a transition layer |29l p. 271], [161 sec. 1.13] of finite thickness 
at an interface between free space and the macroscopic dipolar continuum [iniiiiiiSQ]. Unlike an ideal 
dipolar continuum, all the components of the E and H fields in a macroscopic dipolar continuum 
without the interface are not generally equal to the fields close to an interface between free space and 
the macroscopic continuum. However, the tangential E and H fields with and without the interface 
are nearly continuous across a transition layer of thickness <5 containing the interface. Fortunately, 
under the conditions that the discrete dipolar material behaves as a continuum, namely <C 1 

and <C 1 (see Section (3^ . analytical and numerical results with discrete dipolar arrays indicate 

that the thickness 5 of the transition layer is on the order of a few average separation distances d of 
the dipoles |3n[l31jPI This means that the infinitesimal thickness of the free-space shell about S that 
was chosen to derive the results for the Poynting vectors in Section 14.11 for the ideal continuum can be 
made equal to the thickness 5 of the transition layer without appreciably changing the average fields 
within the volume V. Consequently, the equation (1551) holds to a good approximation for S at the 
center of the 5-thick free-space shell that removes a discrete number of dipoles (see Fig. ED, and the 
total instantaneous power flowing across the closed surface S in free space just outside the volume V 
containing a discrete number of dipoles is given to a good approximation by (|56ll . that is 



fi • [E(r, t) X H(r, t)]dS 


V 


dD(r,t) 




dV. 


(58) 


The better the discrete distribution of dipoles approximates a continuum, the thinner the transition 
layer becomes (5 —?■ 0 for an ideal continuum) and the more accurate is the approximation in (|58l) for 
the macroscopic continuum. 

The most important aspect of (|56p for the ideal continuum and in ()58p for the macroscopic 
continuum is that both the surface and volume integrals in (I56p and (|58p are equal to the instantaneous 
power flow P(t) through the closed surface S into the volume V, exactly for the ideal continuum and 
approximately for the macroscopic continuum, respectively. Although free-space shells are invoked to 
prove (l56P and ([581) . the fields in (l56P and (fSHP are the continuum fields and thus (l56P holds exactly 
within the ideal continuum, and (15811 holds approximately in the macroscopic continuum since n- (E x H) 
is continuous across a free-space/dipolar-material interface under the conditions specified above and 
satisfied by most dipolar metamaterial continua as well as natural dipolar material continua. 

® Computations in front of planar arrays of point dipole sources separated by a distance d show that the fractional variations 

in the fields caused by the localization of the point dipoles are no greater than about 1% for a distance z = d\n front of the plane of the 
array; and the fractional variations rapidly decrease with z > d. In other words, fi • (e x h) at the center of a free-space shell just two 
lattice layers thick {5 = 2d) is equal to fi • (E x H) to an accuracy better than a few percent with the accuracy rapidly improving for 
6 > 2d (and the integrations over S of these two Poynting vectors agree to appreciably greater accuracy). Also, extensive calculations 
with three-dimensional dipole arrays show that they are well approximated by continua as the lattice spacing d becomes less than 
about a tenth of a wavelength m- 
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Figure 2. Volume V (or AV) with its surface S (or AS”) centered in a free-space shell of thickness <5 
enclosing a large number of discrete electric and magnetic dipoles. 


5. ENERGY RELATIONS FOR MACROSCOPIC DIPOLAR CONTINUA 


If the volume V in (|58p is chosen to be a macroscopic volume AV that is electrically small but large 
enough to contain a great number of discrete dipoles, then ([58]l becomes 


AP{t) J h- [E(r, t) X H(r, t)]dS = j 

AS AF 


9D(r, t) 

m 


E(r,t) + 


clB(r, t) 


H(r,t) 


dV. 


(59) 


We can also express the power AP{t) in terms of the microscopic fields and sources within AV. This 
can be done by reinstating the surface AS of AV in the middle of a free-space shell of transition 
layer thickness <5 that removes a discrete number of dipoles, as explained in Section 14.21 where the 
volume AV contains a large number of discrete electric and magnetic dipoles. Since AS is now in free 
space a distance 5/2 away from the nearest dipole, the average macroscopic E and H fields there are 
approximately equal to the microscopic e and h. = h/fiQ fields, as explained in Section 01 Thus, from 
(|59|) and ifTSl) . the power AP{t) can be written as 


AP{t) 


-J h-[e{r,t) xh{r,t)/no]dS 
AS 

j(r,t)-e(v,t)dV+~ I 


AV 


AV 


AV 


Al/ 


eo|e(r,t)p + —|b(r,t)p 
ho 


dV 


8D(M) 


dt 


dP(r,t) dP(r,t) 

^ ^ •E(r,t)-^.M(r,t) 


dt 


dt 


dV 


dV 


1 d r 
'^2di J 

AV 


eo|E(r,t)|2 + —|B(r,t)|2 
ho 


dV. 


(60) 
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This equation reveals that 

I 'dP{r,t) 

AV 


dt 




dV 


i j j(r,t) • e(r,t)dy + J |eo [|e(r, - |E(r, t)|2] + ^ [|b(r, t)|2 - |B(r, t)|2] | dF . (61) 

Ay Ay 


The last integral in (I61j) can be simplified using the identities 

j \e-B\‘^dV= y'(|e|2 + |E|2-2e-E)dy = j {\e\‘^ - \B\‘^)dV 


AV 


AV 


AV 


j \h-B\‘^dV= ydbl^ + lBp-2b-B)dE= ^ (|b|2 - |B|2)dE 
Ay Ay Ay 


(62a) 


(62b) 


which hold because E and B within AV are practically uniform over AV and these primary macroscopic 
fields are defined simply in terms of the microscopic fields as 


E(r eAV,t) = ^ I e(r, t)dV (63a) 

Ay 

B{rGAV,t) = ^ I h{r, t)dV. (63b) 

Ay 

Moreover, the microscopic fields within AV from the dipoles outside of AV are practically equal to the 
macroscopic fields from the dipoles outside of AV. That is, both the microscopic and macroscopic fields 
within AV produced by the dipoles outside AV are approximately equal to the cavity fields in a cavity 
AV of the corresponding ideal continuum and thus (e°“* — E°“*) ~ 0 and (b°“* — B°^*) ~ 0, where the 
superscripts “out” denote fields within AV from dipoles out side AV. Therefore, 


|e-E|2« (64a) 

|b-B|2 « |b“'^-B“'^|2 (64b) 

where the superscripts “ins” denote the fields within AV produced by the dipoles inside AV. By 
definition 

E“®(r e AE,t) = ^ y e”'^(r,t)dE (65a) 

Ay 

B^ir eAV,t) = ^ J W^{r,t)dV. (65b) 

Ay 

Combining (j64p with (I62p . inserting the result into (1611) . then performing the volume integration on the 
left-hand side of (j6ip by using the fact that the macroscopic sources and fields are nearly uniform over 
the electrically small macroscopic volume AV, one obtains 


dP ^ dB 

dt dt 


AE 


j • edV + 


1 d 
2AV di 


AV 


AV 


eole*'"" - E“T + —|b“" - B 
do 


ins 1 2 


dV. 


( 66 ) 
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Integration of this last equation from time to to time t gives 

t 

{ T» V ~\ 

dt' 


to 


9P(r, t') , (9B(r,t') , 

^ E(r,0-W.M(r,0 


dt' 


AI/ 


dt' 


j • edVdt' + 


2AV 


to Ay 


Ay 


eo|e“" - + —|b“" - 


1 

t^o' 


11 


dV 


- io 


AF 


j • edVdt' 


to Ay 


+ 


2AF 


Ay 


eo|e“'^(r,t) -E“^(r,t)|2 + —|b“^(r, t) - B‘°^(r, t)|2 - eo|e‘°"(r, to)P - —|b“^(r,to)P 

ho ho 


(67) 


dV 


because E™®(r, to) = 0 and B‘°®(r, to) = 0 . 

This result shows that the change in macroscopic polarization energy in AF determined by the left- 
hand side of ()67p is not just equal to the work done by the microscopic electric field on the microscopic 
current density (the mechanical energy represented by the double integral on the right-hand side of 
()67p ) because of the additional change in field energy on the right-hand sides of ()67p . We will show 
next that this additional field energy produced by the dipoles inside any one macroscopic volume AF 
is conhned to that AF and thus behaves as an isolated capacitive and inductive energy of the system 
of dipoles inside AF. 

Depending on the properties of the discrete dipoles, the last volume integral in (j67p may or may not 
be greater than or equal to zero. To better understand this, it helps to give a physical interpretation of 
the integral of the “inside” fields in (1671) . The helds (e™®, b”*^) are the microscopic fields produced by the 
discrete dipoles located within the electrically small macroscopic volume AF, whereas (E'^'^jB™®) are 
the fields produced by the corresponding continuum polarizations (P, M) replacing the discrete dipole 
moments within AF. (For a spherical AF, E‘“® « —P/(3eo) and B”*^ « 2/roM/3.) Therefore the fields 
and E”*^ both satisfy the quasi-electrostatic Maxwell equations with discrete electric dipole moments 
and continuous electric polarization in AF, respectively. Likewise, the helds b™® and B“® both satisfy 
the quasi-magnetostatic Maxwell equations with discrete magnetic dipole moments and continuous 
magnetic polarization (magnetization) in AF, respectively. For observation points outside AF, the two 
pairs of helds are approximately equal, that is, e“®(r, t) « E‘“®(r, t) and b™®(r, t) « B“®(r, t) for r ^ AF, 
whereas for observation points within AF, the two pairs of helds can differ greatly. Consequently 



Ay 


eo|e‘-(r,t)-E'-(r,t)|2 + -|b‘-(r,t) 

fo 


eo|e“^(r,t) - E”'^(r,t)|2 + J-|b“^(r,t) 

fo 


B“‘^(r,t)|2 

B“^(r,t)|2 


dV 


dV 


( 68 ) 


where Fx) denotes all space. Thus, the electric-held integral over AF can be viewed as local capacitive 
energy of the electric dipoles within AF and the magnetic-held integral over AF can be viewed as local 
inductive energy of the magnetic dipoles within AF, at every instant of time t: 

AWc{r,t) = l J ^0 [|e“^(r,t) - E“^(r,t)| 2 ] dV (69a) 

Ay 

AWL{r,t) = J / - [|b-^(r,t) - Bi-(r,t)|2] dV (69b) 

2 J ^0 

Ay 
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where the r in AWc(r, t) and AVh;.(r, t) indicates the position of the macroscopic volume AV. The 
subtraction of and from e™® and b™®, respectively, in (I68p removes the fields contributed by 
the truncation of the volume AV at the inner surface of the free-space shell shown in Fig. [2j These 
difference fields produced by the dipole moments inside AV do not extend across AS. The energy in 
these inside fields plus the mechanical energy supplied by the microscopic electric field e acting on the 
microscopic current j in AV through the j -e integral on the right-hand side of (1671) (see Section [2.l.ip is 
equal to the change in macroscopic polarization energy density in AV determined by the first integral 
in (l67p . 

If we divide the electric current j into the electric current je of the electric dipoles and the electric 
current j;, of the magnetic dipoles, such that j = je + jfe, and define electric and magnetic polarization 
energy densities as 


L 

to Al/ AF 

t 


(70a) 


Uh{r,t) = 


where 


with 


AF 


[ ^|b“''(r,t)-B“''(r,t)pdF = [AWj^eir,t) + AWLir,t)] 


to AV 


J jjQ 

AV 


AWje{r,t) = AWj^e{r,t) + AWj^e{r,t) = / / ■ e{r,t')dVdt' 


to AV 
t 

AWj^e{r-,t) = J J je(r,t') • e{r,t')dVdt' 

to AV 
t 

AWj^e{r,t) = j j jb{r,t') ■ e{r,t')dVdt' 

to AV 

then the change in macroscopic polarization energy density in p67p can be rewritten as 


(70b) 

(71a) 

(71b) 

(71c) 


to 




dt' 


dt' 


dt' PS Ue{r, t) - Ue{r, to) + Ub(r, t) - Ub{r, to). (72) 


If the initial energies of formatiorO of the electric and magnetic dipoles do not change or they get 
greater with time (that is, no energy can be extracted from the initial energy of the dipoles), then in 
passive material 

t/e(r, t) - I 7 e(r, to) + I7b(r, t) - Ub{r, to) > 0 (73) 

and we have 


to 


dP(r,t') cIB(r, t') , 

^ E(r,0- ^.M(r,0 


dt' 


dt' 


dt' PS Ue{r, t) - Ueir, to) + Ub{r, t) - Ub{r, to) > 0. (74) 


^^The energy of formation of each dipole is the work needed to differentially and quasi-statically assemble the charges or currents in 
free space from an infinite separation distance. It can be rewritten in terms of the integration in free space over the static electric 
and magnetic field of each electric or magnetic dipole, respectively. The energy of formation of a permanent electric dipole does not 
change if the separation distance of its charges do not change. The energy of a “permanent” magnetic dipole does not change if its 
current does not change. 
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To clearly demonstrate this inequality, chose a spherical volume AV and enclose the dipoles in Al^ with 
a perfect conductor so that the electric and magnetic fields in Al^ will equal and — B“®, 

respectively. Then an external source exerting nonelectromagnetic forces equal to e(r, t) on the dipole 
charge carriers will produce the mechanical energy change equal to AWje and the dipole moments pi 
and rrij in AV that, along with the induced surface charge in the perfect conductor, generate 
and — B™®. The total energy supplied by the external source will be the sum of the mechanical 
energy change in the charge carriers and the energy stored in the fields, that is, the energy density given 
on the right-hand side of (j72p . Since no energy can pass through the perfect conductor enclosing AV, 
the right-hand side of (j72l] must be greater than or equal to zero in passive material where no energy 
can be extracted from any initially existing dipole moments. 

A material that satisfies the inequality in (I74p can be referred to as being “unconditionally passive”. 
An example of a dielectric material satisfying (I74p . so that the material is unconditionally passive, is 
one in which its electric polarization is produced by an initially random distribution of isolated electric 
dipoles that can be modeled by equal and opposite electric charges at the ends of rigid rods with each 
rigid rod connected to a rigid lattice by an untorqued torsion spring that can be linear or nonlinear, 
lossy or lossless. 

In contrast, we will show that “permanent” Amperian magnetic dipoles are not unconditionally 
passive because their initial magnetic dipole moments can be reduced as they align in an external 
magnetic field. To illustrate this, consider a causal paramagnetic material with a constant real 
magnetic susceptibility Xm > 0 over the baseband operational bandwidth |a;| < uq such that 
M = + Xm)] = XfoB/z/o, so that > 0 over the baseband bandwidth. Then, for negligible 

P and the macroscopic fields zero at time Iq, the left-hand side of (j74h . the change in macroscopic 
polarization energy density, becomes 


-—\B{r,t)\‘^dV <0 (75) 

2/xo 

which contradicts (I74p . This contradictory result indicates that the energy change in (I73p and (I74|] 
applied to magnetic materials may be a nonnegative energy only for diamagnetic materials, a conclusion 
proven in Section [5.21 below. 

5.1. Evaluation of the macroscopic polarization energy for electric dipoles 

Nearly all electric polarization can be modeled by either initially randomly oriented molecules with 
permanent electric dipole moments that stay practically fixed in magnitude (like rigid-rod dipoles) as 
they align in the applied field, or by initially zero electric dipole-moment molecules that distort in the 
applied field to produce induced electric dipole moments [8l p. 464], [TTl p. 162]. In the latter case, 
Aiyc'(r, to) = 0 and the electric dipoles have no initial energy of formation that can be reduced. In 
the former case, the energy of formation of the electric dipoles cannot be reduced because the initial 
separation distance between the equal and opposite charges of each dipole cannot decrease to release 
part of their initial energy of formation. (The difference between this former case and the latter case, 
where the change of energy is manifestly greater than or equal to zero because AlTc'(r,to) = 0 and 
both AWcir,t) and AlEjge(r,i) are > 0, is that the baseline energy of formation of the electric dipoles 
in the former case is not zero. However, since this baseline energy of formation cannot be reduced, the 
change of energy in this former case of fixed-magnitude electric dipoles is also greater than or equal to 
zero, as in the latter case.) 

In other words, the change in macroscopic polarization energy for both these models of electric 
dipoles in passive material satisfies the positive semi-definite inequality in (1741) . Felsen and Marcuvitz [71 
sec. 1.5] arrived at a similar conclusion for a linear, lossless, dispersive, anisotropic medium but without 
the proviso requiring unconditional passivity, which we shall see is crucially important when we deal 
with paramagnetic material in the next section. The derivation of the left-hand side of (1741) in terms of 
the microscopic fields on the right-hand side of (1741) has revealed the necessity of assuming unconditional 
passivity for the microscopic dipoles to ensure that this energy is positive semi-definite. 

In summary, assuming either of these two reasonable models for the microscopic electric dipoles in 
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passive material, we find that if the macroscopic magnetization M(r, t) in (17411 is equal to zero, then 

t 

—« Ueir,t)-Ue{r,to) = [AWj^e{r,t) + AWc{r,t) - AVFc(r,to)] > 0. (76) 

to 

It is difficult to conceive of a simple model of electric dipoles that would disobey the nonnegative 
inequality in (I76l) . For example, if the rigid rods separating the equal and opposite electric dipole 
charges were replaced by initially compressed springs, then the potential energy in the springs could be 
released upon being uncompressed by applied fields to produce a negative AWj^e{t) (so the material 
is nominally nonpassive). However, the accompanying increase in the formation energy of the electric 
dipoles caused by the increased separation distance between their equal and opposite charges is greater 
than the energy released by the springs and, thus, the associated change in the capacitive energy, 
AVFc(r,t) — AVFc(r,to)) would add to AWj^e(t) to keep ([76]l nonnegative. Conversely, if the springs 
were compressed further than their initial compression, the decrease in formation energy of the electric 
dipoles (producing a decrease in AWc'(r, t) — AWc'(r, to)) would be more than cancelled by the increase 
in potential energy of the springs (increase in AWj^e{t)) to, again, keep (I76]l nonnegative. It should be 
noted that even for passive material (that is, AWj^eit) > 0) comprised of permanent dipoles with fixed 
magnitudes (such as the rigid-rod model of electric dipoles), Aiyc'(r,t) — AIFc(r,to) is not necessarily 
greater than or equal to zero even though Ue{r,t) — Ue{r,to) > 0. 

5.2. Evaluation of the macroscopic polarization energy for magnetic dipoles 

Because magnetic charge does not exist, magnetic dipoles are created by molecules with circulating 
electric-charge currents that can be modeled by perfectly electrically conducting (PEC)E3 wire loops. 
If the wire loops carry no permanent current so that all the circulating current and magnetic dipole 
moments are induced by the applied helds, then the distribution of molecules forms a diamagnetic 
macroscopic continuum whose low-frequency magnetic susceptibility is less than zeroJ3 If the wire 
loops carry “permanent” currents with magnetic dipole moments that align in an external magnetic held 
and dominate any induced diamagnetization, then the distribution of molecules forms a paramagnetic 
(which includes ferro(i) magnetic and antiferromagnetic [HI chs. 11 and 12]) macroscopic continuum 
whose magnetic susceptibility is greater than zero at low frequencies. 


5.2.1. Macroscopic polarization energy for diamagnetism 

In the case of diamagnetism, the initial value of the currents and magnetic dipole moments of the wire 
loops are zero and AIFi(r,to) = 0. Therefore, passive diamagnetic-material continua that can also 
contain electric polarization produced by electric dipoles with hxed-magnitude or initially zero electric 
dipole moments exhibit the unconditional passivity inequality in (I74p . specihcally 


/ 


'dF{r,t') 

m' 


E(r,t') 


dB{r,t') 


M(r,t') 


dt’ >0. 


(77) 


Since perfect electric conductors (PEC’s) with infinite conductivity and zero internal electric and magnetic fields at all frequencies 
including oj = 0 do not exist naturally except at superconducting temperatures, it would be more accurate to use the term 
“superconductor” rather than “PEC”. Nevertheless, we will use the classical term “PEC” in the ideal sense of a “superconductor”. 
^^The Bohr-Leeuwen theorem, which states that a classical material composed of free charges in thermal equilibrium cannot be 
affected by a magnetic field because the magnetic field does no work on moving charges [231 secs. 24—27], is often used to argue that 
classical electromagnetics cannot describe diamagnetism. This theorem assumes that all the energy at thermal equilibrium resides 
in the kinetic and potential energy of the charge carriers. However, by using a conducting loop model for molecules or metamaterial 
inclusions, an additional inductive magnetic-field energy is introduced such that the Bohr-Leeuwen theorem no longer applies and we 
are able to classically derive macroscopic power energy relations in diamagnetic materials and metamaterials. In fact, W. Weber and, 
to a greater extent, Maxwell, in his Treatise [2] arts. 836—845], explained both diamagnetism and ordinary magnetism by means of the 
perfectly conducting wire loops with no “primitive current” (Maxwell’s term for “permanent current”) in the case of diamagnetism, 
and predominantly “primitive current” in the case of paramagnetism. 
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As a simple example, consider a medium with negligible loss and constant real permittivity e and 
permeability /r over the operational baseband bandwidth |c(j| < ojq. Then P = (e — eo)E and 
M = (l//ro — over the baseband bandwidth, and (f77j) reduces to 

(e - eo)|E|2 - (1/^0 - l//u)|B|2 > 0. (78) 

Since E and B can be chosen equal to zero independently, this inequality reveals that 

e > eo , 0 < fJ-o (79) 

which confirms that the inequality in (I77p applies to diamagnetic continua. 

The inequality in (fTTl) was used in [33], without showing how it was derived, to obtain inequalities 
satisfied by frequency-domain diamagnetic permeability. A derivation of (|77p was first given in |34j . 
Inserting P and M from the constitutive relations (1481) into (I77p recasts this inequality into the form 


t 



to 


dn{r,t') , dB{r,t') 


_1 

ei. 

IV 

3 1 1—‘ 

eo|E(r,t)|2 -h — |B(r,t)|2 

1 

L do \ 


(80) 


which holds for unconditionally passive, spatially nondispersive continua with electric dipolarization 
and diamagnetization. An important feature of (j80p is that the stored macroscopic diamagnetic field 
energy density is |B(r, t)p/(2/ro), which is not generally equal to /io|H(r,t)p/2 (except in free space). 


5.2.2. Macroscopic polarization energy for paramagnetism 


For paramagnetic substances, which include ferro(i)magnetic and antiferromagnetic substances [U chs. 
11 and 12], the molecules or inclusions contain “permanent” magnetic dipole moments that can be 
modeled by randomly oriented, translationally stationary, PEC wire loops carrying initial currents 
/q with magnetic dipole moments mo in the absence of externally applied fields such that when an 
external field is applied, the magnetization produced by the partial alignment of these magnetic dipole 
moments dominates the magnetization produced by the induced diamagnetic currents. The physics 
would appear to be similar to that of fixed-magnitude permanent electric dipoles and, thus, one is 
tempted to apply the argument used in Section 15.11 for permanent electric dipoles to conclude that 
(|77p and (1801) hold also for paramagnetic materials. However, there is an important subtle difference 
with these “permanent” Amperian magnetic dipoles that requires a more involved analysis than for the 
fixed-magnitude permanent electric dipoles. Even though the “permanent” magnetic dipole moments 
can be assumed to dominate over the additional diamagnetic dipole moments induced by applied fields, 
these small additional induced diamagnetic dipole moments can significantly change (and, in particular, 
reduce) the initial energy stored in the magnetic fields of the wire loops, as we shall show. This reduction 
in initial internal energy invalidates the unconditionally passivity argument that was used in Section 
o for fixed-magnitude electric dipoles to show that (I77p holds for all t > t^. In fact, as indicated in 
dZi, the inequality in (I77p does not generally hold for paramagnetic continua. 

Even if we return to (1601) and rewrite the last equation in (|60p as 
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and obtain, instead of (IGTI) . the equation 
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which, when integrated over time gives 

t 


to AF 
t 

-JJ ■e{r,t')dVdt' + ^ j {cq [|e(r, - |E(r, t')P] + /^o [|h(r, t')P - |H(r, t')P] dE (83) 

to Ay Ay 

it is not apparent that the right-hand side of this equation in ([83l) is always greater than or equal to 
zero for paramagnetic macroscopic continua because the initial energy of formation of the Amperian 
magnetic dipoles can be reduced. Also note in (I83p that |hp — |Hp 7^ |h — Hp because for Amperian 
magnetic dipoles f^y \idV = ^^yhdV/ hq = BAV"//io / HAl/. 

To derive a positive semi-definite macroscopic polarization energy for the “permanent” microscopic 
Amperian magnetic dipoles of paramagnetic materials, we begin by evaluating Pje{t) for a microscopic 
(electrically small), perfectly electrically conducting, thin, rigid, wire loop spanning the area A and 
carrying an initial current Jq such that the initial magnetic dipole moment of the loop is mo = loAn, 
where n is the normal to the loop in the direction determined by the right-hand rule applied to the 
circulation of Iq. For the sake of simplicity, assume that the thin wire loop lies in a plane (see Fig. 
H, although the derivation goes through for electrically small wire loops of any shape characterized 
by a self-inductance. The loop can be considered bound to a rigid lattice by lossless or lossy torsion 
springs that allow the loop to rotate about axes through its center rg. With the wire loop illuminated 
by external electric and magnetic fields, eext(r,t) and bext(r,t), we can write 

PjbS) = j ■e{r,t)dV = j j6(r,t) ■ ee^t{r,t)dV + j jfe(r,t) •eind(r,t)dE (84) 

y y y 

where eind(rT) is the induced electric field, that is, the electric field generated by the current jb(r,t) in 
the loop. Note that because the loop can rotate, the total electric field e = Oext + eind is not equal to 
zero in the PEC wire. Rather it is the Lorentzian force field that is zero in the rotating PEC wire, that 
is 

e(r, t)-|-v(r, t) X b(r, t) = 0 (85) 

where v(r, t) is the velocity of the point r in the wire at the time t and the total magnetic field 
can be written as the sum of the externally applied magnetic field and the induced magnetic field, 
b = bgxt + bind- The induced magnetic field bind(r, t) is produced by the current in the loop, including 
the initial current Iq. 

For an electrically small loop carrying a current I{t) at the time t, the integral in (I84|) involving 
Oext can be evaluated quasi-statically as 



dP{r,t' 


• E(r,t') + fio 




•H(r,0 


dVdt' 


J ■ eext{r,t)dV = I{t) j) eext{r,t) ■ dc =-I{t) J —^{r,t)-ndS 

V c(t) S(t) 

w,s<9bext(ro,i) f^,c, Tf^\A- ^bgxtfroT) ... 9bext(roT) 

= -Hi) -- ndS = -I{t)An - — -= -m(t)- — - (86) 

s 

where C{t) and S{t) are the rotating curve and planar surface, respectively, defining the thin wire loop 
at the time t, and c is the unit vector along the direction of the wire. The result in (I86|) expresses the 
power supplied by the externally applied voltage to the current in the PEC wire loop in terms of the 
magnetic dipole moment and the time derivative of the external magnetic field at each instant of time 
t. 

A relatively easy way to evaluate the integral in (1841) involving Ojnd for an electrically small loop is 
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Figure 3. Perfectly electrically conducting (PEC), electrically small, rotating, thin, rigid wire loop 
carrying an initial current Iq and subject to an external field so that the total microscopic current 
distribution is j;,, which produces a magnetic dipole moment m. The total microscopic electric and 
magnetic fields, e and b, are equal to the sums of the external fields and the induced fields produced 
by jft. 


to extend the integration from V to all space 14o, that is 



eind(r, t)dV = / jb(r,t) ■ eind(r, t)dV. 




(87) 


Then, insert j;, from Ampere’s law (Maxwell’s second equation), jf, = V x bind//^o “ eo^Sind/^t, to get 


/ jf)(r,t) • eind(r, t)dV = / 

Voo Voc 


f±Vxbi,d-eo^) -einddV 

V//0 Ot J 


( 88 ) 


which, with the help of the equations, V • (ejnd x bind) = (V x e^nd) • bind — (V x bind) • ^ind and 
V X Oind = —dh\^d/dt, converts to 


^ ^^|bindP + eoleindP^ dV = j ^^|bindP + eo|eindP^ dV 

Voo Voo Voo 

(89) 

after noting that the integral over Soo of the quasi-static Poynting vector is zero. The energy in the 
quasi-electrostatic field of the inductor is negligible compared to the energy in its quasi-magnetostatic 
field and, thus, ([891) can be rewritten as 


jb(r,t) •eind(r,t)dP = / 


j jb(r,t) • eind(r,t)dP 
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dt 



bind(r,i)PdE 
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(90) 


where L is the inductance of the PEC wire loop |35l ch. 8]. 

With the results in (l90]l and (l86l) inserted into (l8^ . we find that the total power supplied by the 
fields to the current loop is given by 



jfe(r,f) • e{r,t)dV 
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(91) 
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This equation (1^ says that the power Pj^eit) supplied by the electric field to the current carrying PEC 
wire loop is equal to the total power supplied by the external fields minus the time rate of change of the 
energy stored in the quasi-magnetostatic field of the current loop. As explained in Section 12.1.11 the 
power Pj^e(t) in (I9ip is also equal to the time rate of change of the mechanical energy, that is, the sum of 
(i) the kinetic energy of the charge carriers and the kinetic energy of the rotating PEC conductors, (ii) 
the kinetic and potential energy of the torsion springs attaching the rotating loop to the rigid lattice, 
and (iii) any heat loss in the torsion springs. 

The inductance power on the right-hand side of (|9ip can also be expressed in terms of the magnetic 
dipole moment m and the externally applied magnetic field bext by invoking the circuit-theory result 
(from the quasi-static Maxwell equations) that the electromotive force applied around the electrically 
small rotating PEC wire-loop circuit by the externally applied fields equals the product of the inductance 
and the time derivative of the current, that is 


[eext(r,t) +v{r,t) x bext(r,t)] ■ dc = L 

C{t) 


(U 

dt 


(92) 


where v(r,t) is the velocity of the thin rotating loop (defined by the closed curve C{t)) at each point r 
of C{t) at the time t. With the aid of Faraday’s law for a moving curve [361 pp. 39-40], namely 


/ d f d 
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so that dSH) becomes 


Pjbeit) = / j6(r,i) ■e{r,t)dV = bext(roA) ' 


dt 


( dm ^ dm 


I dt dt 
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(96) 


in which we have used the identity Adl/dt = dm/dt. Note that Pj^eit) is zero if the wire is not 
rotating (because then hdm/dt = d{mh)/dt = dm/dt), a result that checks with e = 0 in a stationary 
nonrotating PEC wire loop. Since we are assuming that we are dealing with paramagnetic materials 
in which the induced diamagnetic moments are negligible compared with the aligning (rotating) 
“permanent” magnetic dipole moments, it follows that \hdm/dt\ = \dm/dt\ <C jdm/dtj « mojdn/dtj 
and (IM|) reduces to simply 

Pjbe{t) = j jb{r,t) • e{r,t)dV « mobext(roA) ' ~ bext(roA) ' ■ (9'^) 

y 

Although the magnitude of the induced magnetic dipole moment is negligible compared to the magnitude 
of the total magnetic dipole moment as the dipole aligns in an external magnetic field, the change in 
energy produced by the induced magnetic dipole moment is not negligible. 

The result i n (1971) . which was also obtained in [121 eq. (2.164)] for “permanent” current enforced 
on a PEC spheram rather than for the “permanent” current on a PEC wire loop, is quite revealing. It 


Although not stated in na, the electrically small PEC sphere with its static dipole moment has to be rotating. 
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says that the power supplied by the electric field applied to an inclusion or molecule with a “permanent” 
Amperian magnetic dipole moment (modeled by the PEC wire loop carrying the initial current Iq with 
magnetic dipole moment mo ) that predominates over any induced diamagnetic moment (as is apparently 
the case for most if not all known natural paramagnetie materials) is the same as the power supplied by 
the magnetic field to a hypothetieal fixed-magnitude ) magnetic dipole moment formed by separated 
equal and opposite magnetic charge [T2l eq. (2.167)]o 

Under the same condition that the rotating “permanent” magnetic dipole moments dominate the 
induced magnetic dipole moments, we have from ([^511 and pn|) that 


j3b{r,t) ■ eind(r,t)dU 

V 


A 

dt 



— [bext(ro,t) • m] 
dt 


(99) 


which can also be obtained by subtracting (f86|) from dW)) . This power supplied by the internal electric 
field to the current in the PEC wire loop is the microscopic “hidden power” corresponding to the 
time rate of change of hidden momentum m pp. 214-216, 244], [UESIEH] for Amperian magnetic 
dipoles in the electromagnetic force-momentum equation that is analogous to the power-energy equation 
(Poynting’s theorem). Specifically, the time rate of change of hidden momentum corresponding to ()99p 
is equal to —{l/c^)d{ui x eext)/dt [121 eq. (2.161)]. Here, our detailed microscopic derivation has shown 
that the microscopic “hidden power” is drawn from the reservoir of inductive energy in the initial 
microscopic Amperian magnetic dipole moment. 

Referring to natural materials or metamaterials whose magnetism is dominated by the alignment of 
“permanent” Amperian magnetic dipoles as paramagnetic materials, we have shown that paramagnetic 
materials behave energetically as if their magnetic dipole moments were produced by unconditionally 
passive, hxed-magnitude magnetic-charge dipoles (for example, equal and opposite magnetic charges at 
the ends of rigid rods). Although energy can be extracted from the initial internal inductive energy of 
the wire loops (Amperian magnetic dipoles) to invalidate the unconditional passivity used to derive (1771) . 
the external and internal power on the right-hand side of (|91l) combine to form a total power in ()97p 
supplied by the electric field to the Amperian magnetic dipole that is equivalent to the power supplied by 
the magnetic field to an unconditionally passive, magnetic-charge magnetic dipole with the same initial 
magnetic dipole moment. Consequently, a convenient way to determine positive semi-definite energy 
relations for macroscopic paramagnetic dipolar continua that correspond to the energy relations in (I77h 
and (|80p for diamagnetic dipolar continua is to repeat the analysis in Sections [2l through 15.11 beginning 
with the Maxwell microscopic equations that include microscopic magnetic charge, namely [161 p. 464] 


/ \ i9hc(r,f) . , , 

Vxec{r,t) + no = jm(r,t) 

(100a) 

/ X dec(r,t) . , , 

Vxhc(r,t) eo =je(r,f) 

(100b) 

V • hc(r,i) = -^Qm{r,t) 
ho 

(100c) 

eoV • ec(r,f) = £»e(r,f) 

(lOOd) 


where gm{T^,t) and jm(r, t) are the magnetic charge and current densities, and ec(r,t) and hc(r, t) are 
the primary fields in the magnetic-charge magnetic dipole system (with electric-charge electric dipoles). 
It is assumed that any electric and magnetic dipole moments are produced by electric and magnetic 


particular, we have 

j jb(r,t) ■ e{r,t)dV = J jrn{r,t) ■ hc{r,t)dV (98) 

V V 

where }rn is the magnetic current of the rotating fixed-magnitude magnetic-charge magnetic dipole and he is the total (external 
plus induced) magnetic field. (For fixed-magnitude magnetic-charge magnetic dipoles, fyjm • = 0-) Since the electric and 

magnetic fields well outside electrically small dipoles are the same whether they are produced by electric or magnetic charge-current, 
Poynting’s theorem applied to a free-space surface enclosing but well away from either an Amperian or magnetic-charge magnetic dipole 
has the same h- (ex h) = h- (gc x he) and thus the same free-space stored energy fy{eo\e\^-\-\h\^ / fio)dV/2 = fy (eo|ecp+/io |hcP)dV^/2 
because the energy in ll98l l supplied to each of the dipoles is also the same. 
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charge separation (^e and Qm), respectively, and that the electric and magnetic currents, je and jm, 
do not produce Amperian magnetic dipoles or analogous circulating-magnetic-current electric dipoles, 
respectively. Because the power and energy derivations for the paramagnetic dipoles proceed parallel 
to those already given in detail for the electric dipoles, only the main results will be given. 

Instead of (jlSp . we have for the microscopic Poynting theorem that includes microscopic magnetic 
charge-current 

P{t) = - J h - [ec(r, t) X hc(r, t)]dS 
s 

= -heir, t)]dV + ~ J [eo|ec(r, -h/io|hc(r, t)P] dP. (101) 

V V 

Instead of ([20]), we have 

t 

Wjemit) = j j[3eir,t') ■ ecir,t')+]rnir,t') ■h.c{r,t')]dVdt'>0 (102) 

to V 

in a passive material with applied fields zero until after the time to- 

One can replace he with be //tq in (|100p ^ (|102l) because these equations hold in free space. However, 
to derive the conventional form of Maxwell’s macroscopic equations in dipolar continua from the 
microscopic equations in ()100p . he must be retained in (llOOp because now the magnetic dipole moments 
and magnetization are defined in terms of magnetic charge rather than circulating electric current. In 
particular, the equation (|35bp is replaced by 

M(r,t) = — 

AV 

which leads to the conventional form of Maxwell’s macroscopic equations in (14611 if and only if one 
begins with he in (llOOp and not bc//ro- 

Since the power in (llOip is still given by the usual integral of the macroscopic Poynting vector in 
a spatially nondispersive continua, the analysis in Section HI again leads to (l56P and ([581) . However, the 
analysis in Section [5] that led to ([72p , ([77P and (IHOl) for diamagnetic material now leads to the following 
three analogous macroscopic inequalities, (llOdp . ()107p and ()108|1 . for paramagnetic material 

t 


to 

in which 


(9P(r, t) 
dt' 


E(r,t') -F/ro 


5M(r, t) 
dt' 


H(r,t') 


dt' K, [/e(r, t) - Ue{r, to) + Uh{r, t) - Uh{r, to) (104) 


Ue{r,t) = — I I je-ecdPdt'-h / eo|e“"(r, t) - E“"(r, t)pdP 


2AI/ 


to AF 


AV 

1 


[AIP>e.(r,t) + AIPc(r,t)] 


AH 


Uhir,t) = — I I 3n.-h,dVdt' + ^ I /ro|hr(r,t)-H-(r,t)|2dP 


to AV 


AV 

1 


[AWj^hSrp) + ^WLk{r,t)] 


(105a) 


(105b) 
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where 


AW, t) = AW,,e.(r,t) + AW, 


t 

^h^{r,t) = J j \^e{Y,t') ■ec{T,t')+]rn{'f^,t') ■h.c{T,t')]dVdt' (106a) 

to Ay 


with 


t 


^^Tec(rT) = J J 3eir,t') ■ ecir,t')dVdt' 

to Ay 

(106b) 

't 

= J j 3m{r,t') ■ hc{r,t')dVdt'. 
to Ay 

(106c) 


The averaged macroscopic fields and sources are the same for the Amperian and magnetic-charge 
magnetic dipoles. 

The energy density in ()104p is greater than or equal to zero because the initial energy of formation 
of the rigid-rod, magnetic-charge magnetic dipoles (like the rigid-rod electric-charge electric dipoles) 
cannot be reduced to release energy and thus 


to 


9P(r, t') 

dt' 


E(r,t') -h/to 


(9M(r, t') 
dt' 


H(r,t') 


dt' > 0 


(107) 


or 


t 


to 


dT>{r,t') , 9B(r,t') 

^^•E(r,t) + ^^ 


H(r,t') 


dt' > - [eo|E(r,t)p -h^o|H(r,t)p] 


(108) 


in paramagnetically unconditionally passive (with respect to the rigid-rod, magnetic-charge formulation 
of magnetic dipoles), spatially nondispersive material whose applied fields are zero until after the initial 
time to- The two energy equations in (jl07l) and (jlOSp were first obtained in [39lll0] but without the 
rigor of the foregoing derivation of these macroscopic equations from the microscopic equations. Also, 
in [39l|40] it was not explicitly pointed out that these equations do not necessarily apply to diamagnetic 
continua or to spatially dispersive continua with h • (E x H) discontinuous across a free-space/material 
interface; see Section 01 


6. RECAPITULATION OF ENERGY RELATIONS FOR MACROSCOPIC DIPOLAR 
CONTINUA 


Using the macroscopic Maxwell equations and constitutive relations in (1461) and ()48l) . we can rewrite 
the energy relations for unconditionally passive, spatially nondispersive (more generally, continuous 
n- (E X H) across interfaces), electric-dipole/diamagnetic continua as 


V to 


9P(r,T 

W~ 


+ V' X M(r,t') 


B{r,t')dt'dV 


V to 


5P(r,T) 9B(r,t') ,, 

^ E(r,0-^.M(r,0 


dt' 


dt' 


dt'dV > 0 


(109) 
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or 


r z 

j y n ■ [E(r, (') X H(r, = / J 




5 to 


V to 


> 


dt' 


V 


eo|E(r,i)|2 + —|B(r,i)|2 
MO 


dt'dV 


dV 


( 110 ) 


for all volumes V with macroscopic fields zero until after the time to- For the first equation in 
(jlOOp to hold, the surface S of the volume V in ()in9|) must lie in a thin free-space shell enclosing 
the volume V such that the equivalent electric surface currents (delta functions in V' x M(r,t'}) 
are included in the volume integration; see Section 01 (In (|110l) it is irrelevant whether or not S 
lies in free space because n • (E x H) for the material being considered is continuous across a free- 
space/material interface.) Equation (11091) confirms that the equivalent electric volume-current power 
density, Je'*(r, t) • E(r, t), in Maxwell’s macroscopic equations for electric-dipole/diamagnetic material 
is given by [9P(r, t)/dt -|- V x M(r, t)] • E(r, t), that is, Je'^ = dV/dt -|- V x M. Recall that diamagnetic 
material is defined herein as material with magnetization produced by molecules or inclusions that have 
no permanent magnetic dipole moments, only induced magnetic dipole moments. Equations (llOOp - 
(jllOp hold for dipolar metamaterials whose inclusions have no permanent dipole moments as long as the 
frequency is low enough that these metamaterials behave as spatially nondispersive dipolar continua; 
for example, (|109I) ~ (I110I) apply to spatially nondispersive bulk metamaterial continua made of split 
resonant rings |41j . 

Similarly, for unconditionally passive (with respect to equivalent hypothetical magnetic-charge 
magnetic dipoles), spatially nondispersive (more generally, continuous n • (E x H) across interfaces), 
electric-dipole/paramagnetic (including ferro(i)magnetic and antiferromagneti(0) continua 


i 



V to 


'(9P(r, t') 


• E(r,t') /io 


clM(r, t') 

9? 


H(r,f') 


dt'dV > 0 


( 111 ) 


or 

t t 

h • [E(r, t') X H(r, = J J 

S to V to 

>^j [eo|E(r,t)|2 +;,o|H(r,t)|2] dV (112) 

V 

for all volumes V with macroscopic fields zero until after the time to- In either of these equations, it is 
irrelevant whether or not the surface S of the volume V lies in a thin free-space shell enclosing the volume 
V. Equation (lllip confirms that the equivalent electric-plus-magnetic volume-current power density, 
[Je'^(r, t)-E(r, t)-|-Jm (r, i)‘H(r, t)], in Maxwell’s macroscopic equations for electric-dipole/paramagnetic 
material is given by [clP(r, t)/5t • E(r, t) -|- fiodM.{r, t)/dt ■ H(r, t))], that is, J/d = dP/dt and Jm = 
fiodM/dt. Recall that paramagnetic material is defined herein as material with magnetization produced 
by the alignment of “permanent” magnetic dipole moments of molecules or inclusions that dominate any 
induced diamagnetic magnetization. Metamaterials with inclusions characterized by paramagnetic m 
producing a macroscopic paramagnetic magnetization M = f^y r x (V x m)dV/ (2AR) = f^y mdV/AV 
that dominates over any diamagnetism produced by the inclusion’s conduction and electric polarization 
currents would have macroscopic fields satisfying the inequalities in (|llip - (lll2l) . 

Antiferromagnetic materials are comprised of permanent magnetic dipoles that align antiparallel to each other below the Neel 
temperature to produce a paramagnetism in an applied field with small positive-susceptibility, macroscopically-averaged magnetization 
and fields (even though their initial microscopic fields and energies can be large) [Sj ch. 12]. Thus, the macroscopic fields of antifer¬ 
romagnetic materials satisfy the paramagnetic inequalities in (Illll l- (I112|| . 



dD(r,t' 

dt' 


E(r,0 


dB{r,t' 

dt' 


H(r,t') 


dt'dV 
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The difference between the power density integrands in the paramagnetic and diamagnetic continua 
in (lllljl and (|109j) . respectively, is equal to the macroscopic “hidden power” 

(113a) 

with the associated macroscopic “hidden energy” 

(m • B - ^|M|2) = i - ;,o|H|2^ (113b) 

analogous to the “hidden momentum” between the Abraham and Minkowski formulations of the 
macroscopic electromagnetic force-momentum equation |37l pp. 214-216, 244], [F1 I251I55] . As mentioned 
above, our detailed microscopic derivation has shown that the macroscopic “hidden energy” is drawn 
from the reservoir of inductive energy in the initial paramagnetic microscopic Amperian magnetic dipole 
moments. For fields that become periodic after t = to, the time-average of the hidden energy in (I113bp 
approaches a constant. 

The positive semi-dehnite inequalities of the equations in (|lU9p - (|110p and (|llip - (|112p were derived 
from the microscopic equations and thus they have not required linearity or any specihc constitutive 
equations. The functional dependence of the macroscopic polarizations on the helds can be nonlinear 
and anisotropic (or bianisotropic). Both sets of equations, (I109I) - (I110I) for electric-dipole/diamagnetic 
macroscopic continua and (I111I) - (I112I1 for electric-dipole/paramagnetic macroscopic continua, hold to 
a more accurate approximation the better the macroscopic continua approximate ideal continua, that 
is, the smaller the values of and are than unity. Also, a corollary of both these sets of 

equations is that 

t t 

h-[B{r,t') xIiir,t')]dt'dS = j j 

S to V to 

for both diamagnetic and paramagnetic continua. It follows from ()114l) . (I55p . and (IlSp that 



(9D(r, t) 

9? 


E(r, t') -b 


(9B(r,f) 

dt' 


•H(r,0 


dt'dV>0 (114) 


t t 

J J ^ h{r, t')]dt'dS = — J J [6c(rT0 ^ hc{r,t')]dt'dS 

S to S to 

t 

■e{r,t')dt'dV + ^ J [eo|e(r, t')|^ + |b(r, t')| V/«o] dE (115) 


V to 


V 


[je(r,t') • ec(r,t') +j^(r,t') ■ hc(r,t')jdt'dV + - [eo|ec(r, t')P + dV >0. 


V to 


V 


Thus, both the total microscopic and macroscopic electromagnetic energy entering an initially unexcited, 
passive, spatially nondispersive volume of either diamagnetic or paramagnetic dipolar material is always 
nonnegative provided the surface of the volume lies in free space and does not cut through the dipoles. 
It also follows from (11111) that (|83p is greater than or equal to zero for paramagnetic continua. 

Lastly, we note that because the volume-time-integral inequalities in ()109l) - (|110p . (|llip - (|112p . and 
()114p hold for all V, they hold as well with the volume integrals omitted (as long as all averages are 
done over macroscopic volumes), that is, with only the time integrals at each value of position r. In 
addition, as mentioned after (I2UP . the materials to which (|lU9p - (|llUp and (|llll) - (lll2p apply can contain 
conduction currents if the predominant energy change produced by the conduction current is in the form 
of heat (as is usually the case). In other words, the usual conduction current J(r,t) can be added to 
dP(r,t)/dt in (|109p - (lll0ll and (Illll) - (lll2p . 
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7. APPLICATION OF ENERGY RELATIONS TO BIANISOTROPIC CONTINUA 


In this section, the energy theorems summarized in the previous section are applied to linear, passive, 
spatially nondispersive media to obtain frequency-domain expressions for internal energy densities in 
lossless media and inequalities that the linear constitutive relations must obey in lossless media. The 
most general linear, spatially nondispersivconstitutive relations are those for bianisotropic media 
and are given in the frequency domain as 

Da;(r) = e(r) • E^(r) -h T(r) • H^(r) (116a) 

B^(r) = 7l(r) • H^(r) I7(r) • E^(r) (116b) 

where 7l(r), e(r), and [I7(r), r(r)] are the permeability dyadic, the permittivity dyadic, and the magneto¬ 
electric dyadics, respectively. Like the fields, they are, in general, functions of frequency oj and position 
r within the media. 

First, assume the magnetic polarization of the bianisotropic material is produced by paramagnetic 
dipoles. Then equations (I108|) and (I112p apply, that is 


t 


to 


dT>{r,t') , aB(r,t') 


H(r,t') 


dt' > - [eo|E(r,t)p -F ^o|H(r, t)p] 


(117) 


In |39II4L)| , the E and H in (11171) were given a time dependence that began at a zero value and increased to 
a sinusoidal time dependence in order to prove the following inequalities for paramagnetic bianisotropic 
material that is lossless (/I = e = e"*"*, 17 = r"*"*, with superscript “T” denoting the transpose) in 
a finite frequency window about the frequency cj of interest 

Re{E* •(a;e)'-E^ + H* •(a;7l)'-H^ + E^- [(u;(l7T + r*)]'-H*} > [eojE^p + ^o|H^|"] (118) 


[iujeii)' - eo] > uje'ii/2 > 0 (119a) 

[{uniiY - no] > > 0 (119b) 

where the primes denote differentiation with respect to uj and the double index ll indicates the xx, yy, 
or zz diagonal elements of the dyadic. 

If the magnetic polarization of the bianisotropic material is produced by diamagnetic dipoles, then 
equations (IHOl) and (jllOl) apply, that is 


t 



5D(r,t') dB{r,t') 


H(r,t') 


dt'>- 


eo|E(r,t)p -h —|B(r,t)p 

no 


( 120 ) 


^0 

Using an analysis similar to the one in [3M30] but applied to (|120l) instead of (11171) . we find for a lossless 
frequency window 


Re {e: • (we)' • E,, + H: • (w/l)' • + E^ • [(^(17^ + r*)]' • H^} 


> 


no 


( 121 ) 


~ ^o] ^ > 0 (122a) 

^®From Maxwell’s equations, Ho; is proportional to the spatial derivatives in V x Eo? and E^j is proportional to the spatial derivatives 
in V X Hew- Thus, the bianisotropic constitutive relations in some contexts are considered inherently spatially dispersive. Nonetheless, 
the curl spatial derivatives do not excite surface delta functions in or M^^;. Consequently, n • (E^^; x is continuous across 

interfaces and our derived energy relations can be applied to the fields and polarizations that satisfy the bianisotropic constitutive 
relations in llllfill . 
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with a similar (though not identical) inequality for the diagonal elements of the permeability dyadic if 
the orientation of the xyz coordinate system is chosen at each point in space to make the permeability 
dyadic diagonal (which is always possible because the lossless permeability dyadic is Hermitian), namely 

[{oJim)' - yl/yo] > ujyii/2 > 0 . (122b) 

Note that as w ^ 0, ()119all and ()122ap imply that for both paramagnetic and diamagnetic material 

Qz(w —t 0) — eo > 0 (123) 

(assuming lim^_j,o(we(j) = 0). Also, (jll9bp implies that 

0) — /xq > 0 (124) 

for paramagnetic material (assuming hm^_,.o(cu/u(;) = 0). However, ()122bh implies that 

^ 0)[1 — y,ii{u! ^ 0)/ yo] > 0 (125) 

or, equivalently 

0 < yii{uj ^ 0) < yo (126) 

for diamagnetic material (assuming = 0). 

7.1. Subtleties associated with diamagnetic continua 

The inequalities in ()119ap and (|119bll can also be obtained from the Kramers-Kronig causality relations 
for e and /I, respectively [26l sec. 84], Thus, one can ask why (jll9bp does not hold also for diamagnetic 
material or for an array of diamagnetic inclusions that behaves as a continuum. This question was 
answered in detail in mii\ with the answer summarized in the Conclusion of m as, 

“At high enough frequencies w, every natural material or artificial material (metamaterial) no longer 
behaves as a continuum satisfying the traditional time-harmonic dipolar macroscopic Maxwell equations 
with spatially nondispersive constitutive parameters. Although this departure from a continuum be¬ 
havior at high frequencies can often be ignored with impunity for the electric and para/ferro(i)magnetic 
polarization of materials and metamaterials, we show in the Introduction that it is mathematically im¬ 
possible to characterize a material or metamaterial that is diamagnetic and lossless at low frequencies 
cu by a causal [obeying the Kramers-Kronig relations] spatially nondispersive permeability that satis¬ 
fies continuum passivity conditions and whose value approaches the permeability of free space [/xq] as 
the frequency ui approaches infinity. Moreover, this noncausality in the spatially nondispersive dipolar 
continuum description of diamagnetism is more fundamental than the noncausality, discussed in m, 
introduced by the point dipole approximation for scattering from the inclusions (molecules) and is not 
removed by including higher-order multipole moments in the spatially nondispersive continuum formu¬ 
lation of Maxwell’s equations.” 

It is possible for diamagnetic causality to be restored and for (llI9bp to remain valid for a lossless 
diamagnetic continuum if yo is changed to yff = yuiuj —>■ oo) < /xq such that yu{uj —0) — /x“ > 0. 
However, this possibility is unrealistic for spatially nondispersive continua because the magnetic dipole 
moment of a given molecule or inclusion excited by a free-space incident field generally approaches zero 
as a; ^ oo and thus /x“ = /xq [T01[TT]P^ 

Another subtlety arises with regard to a lossless (in a frequency window), homogeneous, spatially 
nondispersive, bianisotropic continuum. In such a continuum, it can be proven [43] that the group 
velocity Vg(ko) of a source-free wave packet (with propagation vectors k concentrated about the vector 

one postulates a spatially nondispersive linear continuum with a hypothetical causal diamagnetic that approaches = 

oo)I, where causality requires that —>• 0) > /i°°, then the method used by Glasgow et al. [4] for continua with causal, 

reciprocal, passive e and implicitly paramagnetic JI that approaches /iol as a; —>• oo can be applied to the hypothetical causal. 
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ko and a narrow-band frequency spectrum t(j(k) such that the envelope of the wave packet is a virtually 
unchanging function g[r — Vfca;(ko)t]) equals the energy-transport velocity Ve(ko); specifically 

Vg(ko) = Vfca;(ko) = = Ve(ko) (129) 

where 

S(ko) = ^Re[Eo X H^] (130a) 

and 

t/(ko) = ^Re{E:-(a;e)'-Eo + H:-(a;7l)'-Ho + Eo- [(^(17^ + r*)]'• H^} . (130b) 

The subscripts “o” on Eq and Hq refer to these fields evaluated at the frequency a;(ko). For paramagnetic 
continua, the inequality in (jllSp shows from (|129l) that 


|Vg(ko)| = |Ve(ko)| < 2 


|Re[EoXH*]| 
eo|Eop + Mo|Hc 


(131) 


Since 


|Re[Eo X H^ll < |Eo||Ho| = -[eolEo^ + " (v^|Eo| - 

< |[eo|Eo|2 + ^o|Ho|2] 


(132) 


where the free-space speed of light is given by c = 1/y^/xoio, this inequality in (jl32p reduces (I13ip to 

|vg(ko)| = |ve(ko)| < c (133) 


confirming that the speed of a lossless wave packet in the paramagnetic continuum cannot be faster 
than the speed of light in free space. 

The subtlety arises by repeating the derivation for a diamagnetic continuum wherein the inequality 
in (jllSp must be replaced by the inequality in (I12ip . so that (I13ip is replaced by 


|Vg(ko)| = |Ve(ko)| < 2 
which can be rewritten from (|132p as 

|Vg(ko)| = |Ve(ko)| < C 


|Re[EoXH*]| 

£o|EoP + I Bop//To 

fo|EoP + ^o|HoP 
€o|Bop + |BoP//io 


reciprocal, passive diamagnetic that approaches fj,°°l as a; ^ oo to obtain the diamagnetic energy inequality 
t t 

r +'\ 1 

dt'dV 


-j j xY-[Y.{r,t') xU{r,t')]dt'dS = j J 


S to 


to 
> 




at' 


at' 


i|[6o|E(r,OP+/.°“|H(r,OP] dV 


from which one finds instead of (I126II for zero loss as a; ^ 0 

^ 0) > 


(134) 


(135) 


(127) 


(128) 

a result that merely repeats the requirement for a causal diamagnetic-continuum permeability. Moreover, neither natural materials nor 
metamaterials are characterized by spatially nondispersive constitutive parameters above frequencies for which the average electrical 
separation distance of the molecules or inclusions becomes a significant portion of a wavelength and, thus, one does not generally have 
a causal expression for diamagnetic material that reveals a reliable analytic high-frequency continuation (in particular, a fj,°° < /^o) 
that can be used in ilmf — (I128II . The inequalities in ilTMf — for electric-dipolar-diamagnetic macroscopic continua and cnj- 
ilTTa for electric-dipolar-paramagnetic macroscopic continua are derived without requiring linearity or reference to any particular 
constitutive relations or parameters. 
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Considering, for example, a lossless isotropic diamagnetic continuum at low frequencies such that 
Bo = /uHo, 0 < fi < Ho, we see that the right-hand side of (|135p can be greater than the free- 
space speed of light c. This result does not actually imply that the speed of the wave packet can be 
greater than the free-space speed of light c in diamagnetic materials, but it is disappointing that the 
right-hand side of (I135p does not come out to be equal to c as it does for paramagnetic materials. What 
it does imply, however, is that the equality in (11291) between the group and energy-transport velocities 
in lossless, linear, spatially nondispersive continua does not provide a robust method for proving that 
the magnitude of these velocities are less than or equal to the free-space speed of light. 

The fundamental principle that can be invoked to obtain this result is the Einstein mass-energy 
relation. Since a fixed amount of energy (call it Wo) in the source-free wave packet of width i travels a 
distance L S> ^ in the lossless continua without significantly changing shape and with velocity v equal 
to the group and energy-transport velocities (v = = Ve), the effective moving and rest masses, m 

and niQ, of this wave packet of energy Wq are found from the Einstein mass-energy relation as 


Wo = mc^ 


moc^ 

x/l-|v|Vc2 


(136) 


which implies 

|v| = |Vg| = |Ve| < c (137) 

for diamagnetic as well as paramagnetic continua. In other words, if a passive medium can propagate 
a fixed-shape bundle of a fixed amount of energy, then the speed of propagation of this bundle cannot 
be greater than the free-space speed of light. A lossless bianisotropic continuum supports just such 
a fixed-energy, distortionless wave packet. Whereas in a lossy, highly dispersive continuum, both the 
amount of energy and the local shape of the wave packet envelope can change and this distortion can 
make parts of the envelope (but not the leading edge of the wave packet) travel faster than the free-space 
speed of light [H] . 

One of the simplest applications of the inequality in (|137l) is to a lossless isotropic magnetodielectric 
material at low frequencies where its e and /r are virtually independent of frequency. Then the group, 
energy-transport, and phase velocities of a low-frequency pulse in the material are approximately the 
same and equal to Ij Thus, (11371) reveals that 


eo “ At 


(138) 


For paramagnetic material, /x > /tq and the inequality in (11381) gives less information than the one in 
(|123p and thus it can be ignored. However, for diamagnetic material, (jl38p can be expressed as 


e > (139) 

where = /xo/l^ is greater than 1 and thus (I139p replaces (|123p . Wood and Pendry [l5] and Sohl et 
al. |46] have also obtained the inequality in ()139p for diamagnetic materials. 


8. CONCLUSION 

With the insight obtained from a relatively simple, yet rigorous proof that averaging the microscopic 
fields of classical discrete dipoles over macroscopic volumes leads to the same Maxwell equations and 
cavity fields obtained for mathematically defined fields in ideal continuous dipolar media, positive 
semi-definite expressions are derived for macroscopic time-domain energy density in passive, spatially 
nondispersive (more generally, continuous n • E x H across interfaces) dipolar continua. The derivation 
proceeds from the underlying microscopic Maxwell equations satisfied by the microscopic fields of the 
electric charge and current that produce the distribution of discrete, bound, dipolar classical models 
of the molecules or inclusions comprising material or metamaterial continua. The bound microscopic 
electric dipoles are assumed to have either zero dipole moments before any external fields are applied 
or randomly oriented, fixed-magnitude dipole moments that can be aligned by the external fields. In 
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both cases, the applied fields can extract no energy from the initial electric dipoles and thus the electric 
dipoles are “unconditionally passive”. 

The microscopic derivation reveals two distinct positive semi-dehnite (nonnegative) macroscopic 
time-domain energy expressions, one that applies to diamagnetic continua and another that applies 
to paramagnetic continua, which includes ferro(i)magnetic and antiferromagnetic materials. The 
diamagnetic dipoles are unconditionally passive because their Amperian magnetic dipole moments are 
zero in the absence of applied fields. The analysis of the paramagnetic continua, which are defined 
in terms of magnetization caused by the alignment of randomly oriented permanent (yet slightly 
changeable in magnitude as they rotate in an external magnetic held) Amperian magnetic dipole 
moments that dominate any induced diamagnetic magnetization, is greatly simplihed by hrst proving 
that the microscopic power equations for rotating “permanent” Amperian magnetic dipoles (which are 
not unconditionally passive because energy can be extracted from their initial magnetic dipole moments) 
reduce effectively to the same power equations obeyed by rotating unconditionally passive magnetic- 
charge magnetic dipoles. 

The difference between the paramagnetic and diamagnetic energy expressions is equal to a “hidden 
energy” that parallels the hidden momentum often attributed to Amperian magnetic dipoles. It is 
noteworthy that the microscopic derivation reveals that this hidden energy is supplied by the change 
of inductive energy in the Amperian magnetic dipole moments as the “permanent” dipoles align in 
an applied held, even though the magnitudes of the induced magnetic dipole moments are negligible 
compared to the magnitudes of the total magnetic dipole moments. 

The macroscopic, positive semi-dehnite, time-domain energy expressions are applied to lossless 
bianisotropic media to determine the inequalities obeyed by the frequency-domain bianisotropic 
constitutive parameters, namely the permittivity, permeability, and magneto-electric dyadics. As one 
would expect, the inequalities obeyed by diamagnetic and paramagnetic permeabilities of bianisotropic 
media are appreciably different. Subtleties associated with the causality as well as the group and 
energy-transport velocities for diamagnetic media are discussed in view of the diamagnetic inequalities. 
In particular, it is shown that general proofs of the group and energy-transport speeds being less than (or 
equal to) the free-space speed of light in lossless media must rely ultimately on the Einstein mass-energy 
relation. 

Lastly, it should be noted that materials or metamaterials have not been considered in which 
their molecules or inclusions consist of a combination of induced diamagnetic dipoles (zero initial 
Amperian magnetic dipole moments) and paramagnetic dipoles (predominant “permanent” magnetic 
dipole moments) such that both the diamagnetism and paramagnetism are comparable. A metamaterial 
array with inclusions containing ferromagnetic material and electric conductors (and/or dielectrics) 
creating diamagnetism would be one such possibility. Then, neither the diamagnetic nor paramagnetic 
energy expression derived in this paper would necessarily apply to this combined magnetization. The 
generalization of positive semi-dehnite macroscopic energy expressions to materials with comparable 
combined paramagnetic and diamagnetic magnetization is a subject for future research. 
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